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ABSTRACT 


This  report  describes  21  Investigations  In  the  field  of  seismic 
discrimination.  These  are  grouped  as  follows:  surface-wave 
studies  (6  contributions),  studies  of  seismic  scaling  and  body- 
wave  magnitude  (4  contributions),  miscellaneous  studies 
(5  contributions),  and  investigations  related  to  data  analysis 
and  computer  systems  (6  contributions). 
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SUMMARY 


1'his  is  the  twenty -seventh  Semiannual  Technical  Summary  report  describing  the  activities 
of  Lincoln  Laboratory,  in  the  field  of  seismic  discrimination.  This  report  covers  the 

period  1 April  to  10  September  1977.  The  objective  of  the  Lincoln  Laboratory  program  is  to 
carry  out  fundamental  research  into  the  seismologies!  problems  associated  with  the  detection, 
location,  and  identification  of  earthquakes  and  nuclear  explosions,  and  with  the  estimation  of 
explosion  yields,  in  order  to  Investigate  these  problems,  we  are  continuously  improving  our 
ability  to  manipulate  and  process  seismic  data  from  a global  network  of  high-quality  digital 
stations  and  arrays. 

Several  investigations  related  to  surface-wave  propagation  are  described.  Polarization 
filtering  has  been  shown  to  be  capable  of  enhancing  the  Rayleigh  waves  received  at  Mashhad 
Seismic  Research  Observatory  (SRO)  from  presumed  explosions  in  Eastern  Kazakh.  In  partic- 
ular, it  is  demonstrated  that  this  technique  will  be  valuable  in  separating  the  surface-wave 
trains  from  interfering  events.  A preliminary  attempt,  to  use  Mashhad  M0  values  in  application 
of  the  Mfl-mb  criterion  to  Eurasian  events  is  quite  successful.  Separation  of  the  explosion  and 
earthquake  populations  is  generally  good,  but  one  presumed  explosion  is  not  Identified  — perhaps 
due  to  unusual  tectonic  strain  release.  We  have  continued  to  improve  our  analysis  of  the  cruBtal 
structure  at  Novaya  Zemlya  by  comparing  the  data  from  two  test  sites  on  the  island  as  observed 
at  ALPA.  A crustal  thickness  of  45  km  appears  to  give  a best  fit  to  the  data.  A comparison  of 
focal  mechanisms  determined  by  a recent  surface-wave  spectral  amplitude  fitting  procedure 
with  those  obtained  from  long-period  (LP)  body  waves  shows  rather  large  discrepancies.  Our 
planned  experiment  to  carry  out  a large-scale  inversion  of  Eurasian  group-velocity  dispersion 
data  is  continuing,  and  requires  data  from  more  SRO  sites  soon  to  be  operational.  Data  from 
a limited  area  in  Southeastern  Asia,  as  observed  at  Mashhad  SRO,  again  document  the  low  group 
velocities  associated  with  the  Tibetan  plateau,  and  indicate  that  these  velocities  persist  to  the 
north  of  the  plateau. 

We  have  continued  our  efforts  to  understand  the  scaling  of  the  seismic  Bource  and  its  effect 

on  determinations  of  m^.  One  study  has  shown  that  observations  of  mb  vs  period  for  the  Rat 

Island  sequence  (1965)  are  consistent  with  a relatively  constant  stress  drop  of  about  25  bars. 

Estimates  of  seismic  moment  appear  to  agree  well  with  other  methods.  Another  investigation 

has  studied  the  effect  of  the  corner  frequency  and  high-frequency  fall-off  in  the  seismic  spectrum 

on  the  dominant  period  as  observed  through  several  short-period  (SP)  seismograph  systems.  It 

is  shown  that  m^,  as  usually  measured,  is  not  a consistent  measure  of  spectral  amplitude  at 

1 Hz  unless  corrections  for  corner  frequency  and  path  attenuation  are  included.  Information 

about  the  seismic  source  can  also  be  obtained  from  the  saturation  of  the  scale  as  observed 

in  frequency-magnitude  data  if  an  earlier -derived  linear  relation  between  log  (frequency)  and 

log  (moment)  is  assumed  to  be  correct.  Observational  data  appear  to  be  consistent  with  a hlgh- 

-2  -3 

frequency  decay  of  the  seismic  spectrum  intermediate  between  u>  and  w . 

Several  other  studies  have  been  completed.  We  have  carried  out  deconvolution  and  max- 
imum entropy  spectral  analysis  of  the  Mashhad  SRO  recording  of  the  Eastern  Kazakh  event  of 
20  March  1976.  This  was  an  unusual  event,  and  is  clearly  multiple.  Our  best  Interpretation 
of  this  record  suggests  that  at  least  one  of  the  sources  was  an  earthquake,  although  further 


v 


data  arc  needed  to  completely  dispose  of  the  possibility  that  it  was  a sequence  of  explosions. 

We  have  begun  a study  of  discrimination  at  near-regional  distances  by  compiling  a set  of  LRSM 
data  from  U.S.  explosions.  The  effects  of  crustal  inhomogeneity  and  the  difficulty  in  measuring 
the  onset  times  of  S phases  are  apparent  in  this  data  set.  In  another  study,  the  application  of 
the  master-event  technique  to  the  teleseismlc  determination  of  focal  depth  is  examined.  In  the 
study  of  an  aftershock  sequence  in  Iran,  it  was  possible  to  determine  relative  epicenters  and 
focal  depths  to  considerable  precision. 

In  the  general  areas  of  data  analysis  and  computer  systems,  several  studies  are  described. 
We  have  completed  a careful  study  of  the  response  characteristics  of  the  SRO  instrument.  We 
note  that  the  anti-alias  filter  has  been  removed  from  the  SP  SRO  system,  and  we  recommend 
that  attempts  be  made  to  design  a new  anti-alias  filter.  We  have  examined  the  detection  thresh- 
olds on  the  existing  SRO  SP  systems  by  evaluating  the  detections  at  each  station.  We  conclude 
that  several  of  these  thresholds  are  set  too  high.  A new  algorithm  for  the  automatic  picking  of 
first  arrivals  is  described  which  appears  to  show  distinct  advantages  over  current  techniques. 

We  are  continuing  to  upgrade  the  capabilities  of  our  PDP-11  computer  system.  Currently, 
the  major  emphasis  of  our  software  program  is  aimed  toward  implementation  of  a multiple-CPU 
operating  system,  the  development  of  a graphics  package  oriented  specifically  for  seismologlcal 
research,  and  the  Improvement  of  our  interface  with  the  datacomputer. 

M.  A.  Chlnnery 
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i.  SKRKAC  K-WAVK  STL  OIKS 

A.  KNILANC'KMKNT  OK  RAYI.KIGII  WAVKS  OK  EXPLOSIONS 

l' SING  POLARIZATION  FILTERING 

Two-dimensional  polarization  filtering*  is  being  used  to  enhance  Rayleigh  waves  from 
explosions.  Here,  we  show  four  examples  of  enhanced  surface  waves  from  presumed  explosions 
in  Kastern  Kazakh  recorded  at  the  Mashhad  SRO  (MAIO). 

The  polarization  scheme  is  based  on  the  particle  motion  ellipse  of  the  two  components  of 
motion,  Z and  R.1  Kor  Rayleigh  waves,  the  Z and  R components  are  90°  out  of  phase  for  all 
frequencies.  The  method  we  use  is  to  shift  the  phase  of  R by  +90°,  thus  making  all  Rayleigh 
components  in  phase  with  the  Z component.  The  motion  then  appears  to  be  highly  polarized  in 
a linear  direction,  such  that  <1^  » d.,,  where  dj  and  d^  are  the  major  and  minor  ellipse  axes  of 
the  particle  motion  in  the  (R, Z)  plane.  The  major  component  along  d,  is  projected  onto  the  R 
and  Z axes  and  scaled  by  a nonlinear  filter  f which  was  previously  defined  by 

f = 1 - e (1-1) 

where  e equals  the  ellipticity  of  the  particle  ellipse  given  by  c = d2/dj.  finally,  R is  shifted 
by -90°  to  restore  its  original  polarity. 

This  filter  passes  waves  with  amplitude  1 when  d2/d^  = 0 and  attenuates  them  gradually  as 
d2  — dj.  Figure  1-1  shows  this  filter  response  as  a function  of  e.  It  smoothly  attenuates  waves 
as  e — 1,  with  the  sharpest  cutoff  occurring  near  e = i. 

We  have  Incorporated  a better  filter  into  the  polarization  program  which  Is  analogous  to  the 
Butterworth  filter  of  order  n used  in  frequency  filtering.  It  has  the  iormula 

g=[l  Me/eo)nfl/Z  (1-2) 

and  is  displayed  in  Fig.  1-1  for  orders  n = 3 and  6 with  a 3-dB  cutoff  ellipticity  eQ  = 0.2.  This 
filter  produces  better  rejection  above  eQ  than  is  possible  using  Eq.  (1-1).  For  e < eQ,  a high- 
order  filter  produces  a flat  response  for  passing  linearly  polarized  waves  with  some  noise 
present. 

This  filter  was  used  In  the  polarization  program  to  pass  surface  waves  from  four  explosions 
in  Kastern  Kazakh  recorded  at  MAIO.  The  results  are  shown  In  Figs.  1-2,  1-3,  and  1-4.  In  each 
of  these  figures  the  two  top  traces,  (a),  are  Z and  R components,  bandpass-filtered  with  a 
phase-free  6-pole  filter  to  pass  0.02  tc  0.08  Hz.  The  next  two  traces,  (b),  in  each  figure  are 
the  results  of  polarizing  the  data  of  the  top  two  traces  using  the  filter  g of  order  3,  which  is 
displayed  as  a function  of  time  In  trace  5.  The  two  traces  which  follow,  (c),  are  the  polarized 
results  using  a 6th-order  filter,  which  is  displayed  in  trace  8. 

In  Fig.  1-2,  the  explosion  '..as  an  mb  of  5.8  with  a very  clear  Rayleigh  wave  on  the  unpolarlzed 
traces.  The  main  effect  of  polarization  is  to  reduce  the  noise  level  somewhat  without  distorth-g 
the  frequency  content  or  amplitude  of  the  surface  wave. 

Figure  I-3ta)  contains  a smaller  Rayleigh  wave,  masked  in  the  coda  of  a larger  earthquake. 
The  earthquake  coda  is  complex,  containing  many  scattered  arrivals  which  are  not  simply  sur- 
face waves  along  the  R direction.  Thus,  polarization  in  Figs.  I-3(b)  and  I-3(c)  clearly  passes 
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tii.  ('Apit  ion  Kayleign  n.ivc  with  littli'  i-mplitude  reduction  while  reducing  the  earthquake  coda 
level.  Some  dispersion  is  also  preserved  in  the  passed  waves. 

tn  1-  ig.  [-*1(3),  two  reported  Kazakh  explosions  occurred  within  5 min.  of  each  other,  both 
surface  waves  being  masked  by  a larger  earthquake  coda.  The  polarization  results  in  Figs.  I— 4(b) 
and  t-4(c)  show  an  increase  in  the  signal/noise  ratio  (SNR)  of  the  surface  waves  for  the  two 
events,  at  the  expense  of  significant  amplitude  attenuation.  The  filter  values  g are  only  about 
~0.2  5 for  these  two  Rayleigh  waves;  however,  the  filter  peaks  indicate  the  presence  of  surface 
waves  which  could  be  correlated  by  an  analyst  to  clearer  arrivals  at  other  stations. 

These  results  show  that  polarization  filtering  can  be  a useful  detection  and  enhancement 
scheme  for  surface  waves  which  docs  not  introduce  severe  signal  distortion  when  the  SNR  is 
high.  The  sensitivity  of  the  mcbiod  also  makes  it  susceptible  to  false  alarms,  which  would  have 
to  be  eliminated  by  a network  correlation  of  all  the  phases. 

C.  W.  Frasier 

R.  A COMPARISON  OF  LOVK-  AND  RAYLEIGH-WAVE  DISCRIMINANTS 

AT  MASHHAD  SRO 

Theoretically,  the  relative  generation  of  Love  and  Rayleigh  waves  is  a powerful  discrimi- 
nant between  earthquakes  and  explosions.  The  fact  that  Love  waves  are  generated  by  explosions 
at  all  Indicates  some  degree  of  tectonic  stress  release  resulting  in  dipole  and  quadrupole  source 
components  in  addition  to  the  monopole  components  generated  by  an  ideal  explosion  source.  A 
convenient  way  of  representing  this  tectonic  stress  release  is  given  by  the  T’  factor,2  being  the 
relative  strength  of  multipolar  to  monopole  source  components.  Values  of  F as  high  as  3 have 
been  observed.2  Spectral  ratios  of  Love -to -Rayleigh  waves  have  been  suggested  as  a discrimi- 
nant, but  have  not  met  with  much  success.2,3 

We  present  here  the  results  of  a preliminary  test  of  the  relative  discrimination  capabilities 

D 

of  Mfl  vs  mb,  where  Mfl  is  calculated  from  the  maximum  amplitudes  of  Rayleigh  (Mg  ) and  Love 
(M,^)  waves  recorded  at  MAIO.  Our  explosion  dataset  is  necessarily  rather  small  since  we 
only  have  data  for  the  first  9 months  of  1976.  During  this  time  period,  there  were  9 presumed 
explosions  in  the  USSR  whose  epicentral  parameters  are  given  in  Table  1-1.  In  addition,  an 
event  very  close  to  the  Eastern  Kazakh  test  site  took  place  on  20  March  1976;  it  has  been  iden- 
tified as  an  earthquake  in  the  PDF.  A special  study  of  this  event  was  made  in  the  previous 
SATS.1 

Our  earthquake  dataset  consists  of  67  earthquakes  within  Eurasia  recorded  l*.  distances  of 
up  to  50°  from  MAIO.  Figure  1-5  shows  the  location  of  the  sources  used.  For  both  the  earth- 
quakes and  the  explosions,  the  horizontal  seismograms  were  rotated  into  the  epicenter  azimuth 
and  M ^ measured  from  the  transverse  component  through 

t:  i.66  log  A + log(-ijr) 

where  A is  the  maximum  peak-to-peak  amplitude  of  Love  waves,  and  T is  the  dominant  period. 
M,?  was  measured  from  the  vertical  component.  T ranged  from  16  to  22  sec  for  both  explo- 
sions and  earthquakes;  these  periods  for  explosions  differ  considerably  from  those  recorded  on 
WWSSN  instruments  for  which  T Is  often  as  small  as  12  sec.  This  difference  is  primarily  due 
to  the  longer  period  response  of  the  SRO  instrumentation. 


TABLE  1-1 

PDE  EPICENTRAL  PARAMETERS  OF  PRESUMED  EXPLOSIONS  IN  USSR, 
JANUARY  - SEPTEMBER  1976 
(Event  10  It  Identified  by  the  PDE  at  on  Eorthquake) 

Event 

Dote 

Time 

Location 

mb 

1 

15  Jonnary 

04:46:57.6 

E.  Kazakh 

5.2 

2 

21  April 

04:47:57.6 

E.  Kozokh 

5.1 

3 

21  April 

05:02:57.4 

E.  Kazakh 

5.3 

4 

19  May 

02:56:57.9 

E.  Kazakh 

5.0 

5 

4 July 

02*56:57.7 

E.  Kazakh 

5.8 

6 

23  July 

02:32:57.9 

E.  Kazakh 

5.1 

7 

29  July 

04:59:57.5 

W.  Kazakh 

5.8 

8 

28  Augutt 

03:56:57.5 

E.  Kazakh 

5.8 

9 

29  September 

02:59:57.4 

Novaya  Zemlya 

5.8 

10 

20  March 

04:03:39.3 

E.  Kazakh 
(Earthquake?) 

5.1 

Figure  1-6  shows  the  Mg1  and  M^  values  obtained;  it  can  be  seen  that  the  line  = M<? 

Is  a very  good  fit  to  the  explosion  population.  All  the  explosions  fall  below  thiB  line,  but  not 

sufficiently  well  to  be  useful  as  a discriminant.  Figures  1-7  and  1-8  are  M -vs-m.  plots  for 
K I 8 o 

Mg  and  Mg  , respectively.  The  values  are  taken  from  the  PDE;  for  nearly  all  the  earth- 
quakes, and  all  the  explosions,  these  are  based  on  10  or  more  observations.  Due  to  large  inter- 
fering events,  w c were  unable  to  measure  Ms  for  events  1 and  4. 

It  can  be  seen  that,  with  one  exception  (event  6),  the  separation  of  the  explosion  and  earth- 

t R 

quake  populations  is  good.  Mg'  vs  appears  to  be  a slightly  better  discriminant  than  Mg  vs 

m^;  excluding  event  6,  the  separation  of  the  two  populations  is  0.5  Ma  units  using  M<?,  and 
0.8  units  using  M^.  The  failure  of  either  to  place  event  6 in  the  explosion  population  is  dis- 
tressing; this  event  must  have  Involved  abnormally  high  tectonic  stress  release.  This  study 
will  be  continued  ns  more  SRO  data  become  available. 

K.  G.  North 


C.  A RAYLEIGH-WAVE  STRUCTURE  FOR  NOVAYA  ZEMLYA: 

RESULTS  FROM  MORE  DATA  ANALYSES 

The  continuing  series  of  large  (presumed)  explosions  at  Novaya  Zemlya  generated  unusually 
clear  Rayleigh  waves.  Recording.?  of  the3e,  made  at  the  Alaskan  Long- Period  Array  (ALFA), 
were  used  to  infer  the  crust  and  upper-mantle  structure  of  Novaya  Zemlya.  For  this  experiment, 
we  had  data  from  five  explosions  which  appeared  to  originate  In  two  source  regions.  The  first 
Is  a no1  them  test  site  locked  near  the  channel  separating  the  southernmost  Island  in  the  chain. 
The  second  is  a southern  test  site  at  the  southern  end  of  that  Island.  A great-circle  path  pass- 
ing through  these  two  teat  sites  also  passes  very  close  to  ALPA.  Thus,  pairs  of  events  In  the 
northern  and  southern  test  site  s,  rcnpi  ct'  vely,  cun  lx*  used  to  determine  the  structure  In  the 
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TAB  UE  1-2 
PDE  EPICENTERS 


iourcc  region,  i.r.,  t:i<  I'ou-nioii  path  between  the  events  across  the  southern  islaitd  of  Novaya 

4 

/.emlvH.  The  approjn  iao-  method  is  the  single-station,  two-event  method  due  to  Alexander. 

Al  I ‘ A center  sen.  • . see- imigrants  from  the  five  events  are  shown  In  Fig.  1-9.  Correspond- 
ing I’DK  information  is  given  in  Table  1-2.  In  order  to  insure  uniform  epicenter  accuracy,  the 

five  events  were  relocated  using  a common  set  of  35  stations.  After  constraining  the  depth  of 

2 

each  event  to  zero,  the  resulting  epicenter  confidence  error  ellipses  were  all  1.64  km  . The 
relocation  procedure  was  necessary  because  the  resulting  phase  velocities  are  Inversely  pro- 
portional to  the  relative  event  separation.  The  average  test-site  separation  was  approximately 
289  km. 

The  method  consists  of  computing  the  urray -averaged  cross-power  spectral -density  func- 
tion. The  phase  of  this  is  the  phase  velocity  dispersion  in  the  source  region  because  the  explo- 
sive source  phase,  instrument  phase,  and  common  path  phase  velocity  dispersion  from  the 
northern  event  to  AI.I’A  all  cancel.  The  amplitude  of  the  cross  power  is  the  relative  cross  power 
in  the  two  events.  The  phase  velocity  dispersion  difference  between  northern  and  southern  events 
can  be  easily  seen  in  Fig.  1-9  where  the  northern  are  events  1,  3,  and  5,  and  the  southern  are 
events  2 and  4.  The  third  event,  23  August  1975,  is  unusual  In  that  it  has  noticeably  less  LP  en- 
ergy than  the  other  two  northern  events.  It  has,  however,  the  same  apparent  phase  velocity 
dispersion. 

Al.PA  is  Ideally  suited  to  observing  events  in  Novaya  Zemlya.  The  propagation  path  crosses 
the  Arctic  Ocean  basin,  conveniently  avoiding  the  North  American  and  Eurasian  continental 
shelves.  Thus,  Interference  due  to  multipathlng  is  minimal.  Little  of  this  effect  1b  apparent  in 
the  seismograms.  The  instrument  response  and  propagation-path  attenuation  were  such  that 
dispersion  could  be  measured  over  three  octaves  of  signal  bandwidth:  from  8 to  64  sec  period. 

The  measured  phase  velocities,  their  estimated  errors,  and  the  phase  velocity  dispersion 
from  the  best-fitting  model  are  shown  in  Fig.  1-10.  The  rras  error  of  the  input  data  was 
0.082  km/sec.  The  phase  velocity  dispersion  from  the  model  fits  the  data  with  an  rms  error  of 
0.045  km/sec.  Its  layer  parameters  are  given  in  Table  1-3  and  plotted  In  Fig.  I— 1 1.  As  can  be 
seen  in  Fig.  1-10,  the  data  show  slight  phase  velocity  minima  at  approximately  the  36,  21,  and 
13  sec  periods.  This  could  be  due  either  to  multipathlng  interference  or  the  presence  of  dipping 
layers  iri  the  source  region.  The  model  does  not  account  for  either  of  these  effects. 

The  model  was  obtained  by  generalized  linear  inversion  of  the  phase  velocity  dispersion 
data.  The  method  has  been  well  reported  In  the  literature;  our  particular  version  is  due  to  Rodl 
et  al_.'>  The  well-known  Gutenberg  Continent  structural  model6  was  used  as  a starting  model. 
Several  runs  of  the  inverse  program  were  made  with  different  crustal  thicknesses.  A 45-km- 
thick  crust  appeared  to  give  the  beet  fit  to  our  data. 

Ray'elgh-wave-dlsperslon  data  are  known  to  delineate  the  shear-wave-velocity  profile  much 
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better  than  either  the  compresslonal  velocity  or  density  profiles.  We  therefore  consider  our 
shear-wave-velocity  profile  to  be  the  major  result  of  this  inversion  experiment;  it  is  shown 
plotted  in  Fig.  1-12  along  with  the  corresponding  shear-wave  velocity  from  the  starting  Gutenberg 
Continent  model.  The  major  differences  consist  of  lower  velocities  in  the  first  two  crustal  lay- 
ers and  the  virtual  elimination  of  the  low-velocity  zone. 

There  is  no  specific  Information  in  the  open  literature  concerning  the  structure  of  Novaya 
Zemlya.  Kosmlnakaya  et  al.  give  a Moho  depth  contour  map  of  the  USSR  Indicating  that  Novaya 
Zemlya  Is  a northern  extension  of  the  Ural  mountain  belt.  Various  geological  maps  of  the  USSR 
support  this.  The  actual  Moho  depth,  or  crustal  thickness,  given  on  the  contour  map  is  between 
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35  and  40  hm.  Our  data  do  not  support  such  a thin  cruet.  In  fact,  the  original  Gutenberg  Con- 
tinent cruet  is  38  km  tiiick.  it  seems  unreasonable  that  even  an  old  mountain  belt  would  only 
be  aa  t'.ick  as  a typical  continent.  Our  crustal  thickness,  45  km,  is  only  7 km  thicker  than 
Gutenberg's  typical  continental  crust. 

The  lack  of  a low-velocity  zone  in  the  shear- wave -velocity  profile  is  indicative  of  an  old, 
stable  structure.  Our  particular  profile  does,  in  fact,  resemble  Ukranian  shield  profiles  pre- 

g 

sented  by  Kosminskaya  et  gL  Although  these  are  appropriate  to  the  continental  USSR,  they  do 
represent  the  crust  in  the  Ural  mountains  region.  To  summarize:  our  phase  velocity  disper- 
sion data  Indicate  a relatively  old,  stable  shear-wave  structure  with  a 45-km-thlck  crust. 

D.  W.  McCowan 
P.  Glovert 
S.  S.  Alexander' 

D.  A COMPARISON  OF  P-WAVE  AND  SURF  ACE-WAVE:  FOCAL-MECHANISM 

SOLUTIONS  FOR  THREE  SINKIANG  EARTHQUAKES 
o 

In  a recent  report.  Sun  and  Shaub  use  a surface-wave  spectral  amplitude  fitting  procedure 
to  estimate  the  focal  mechanism  and  depth  of  19  events  in  the  1974  Sinklang  earthquake  series. 

As  a test  of  this  procedure,  their  solutions  for  three  of  the  largest  events  were  superimposed 
on  lower- hemisphere,  equal-area  plots  of  P-wsve  first  motions.  P-wave  data  were  obtained 
from  LP  vertical  WWSSN  seismograms.  The  eplcentral  data  and  source  parameters  for  these 
three  events  are  given  in  Tables  1-4  and  !-5.  Surface-wave  mechanisms  and  P-wave  first  mo- 
tions are  compared  in  Figs.  1-13.  1-14,  and  1-15. 

As  the  figures  show,  significant  discrepancies  exist  between  surface-wave  solutions  and  the 
first  motions.  Large  numbers  of  compressions  are  seen  in  quadrants  that  should  contain  only 
dilatations  if  the  surface -wave  solutions  were  compatible  with  the  P-wave  data.  Furthermore, 
the  solution  displayed  in  Fig.  1-15  seems  to  contradict  what  Is  known  about  the  tectonics  of  the 
region.  Shirokova,10  Das  and  Filson,11  and  Molnar,  Fitch,  and  Wu12  have  obtained  many  focal - 
mechanism  solutions  for  this  area  and  have  found  that  the  P-axls  commonly  trends  approximately 
N-S.  The  trend  of  the  P-axls  in  Fig.  1-15  is  nearly  E-W. 

Figure  1-13  also  shows  a P-wave  first-motion  solution  for  event  1.  Parameters  for  this 
solution  are: 


Strike 

Dip 

Slip 

Plane  1 

50* 

66SE 

25' 

Plane  2 

310* 

68NE 

155* 

Trend 

Plunge 

P-axls 

359° 

02* 

T-axls 

269* 

34* 

B-axls 

92* 

57  * 

Plane  1 Is  well  defined,  and  plane  2 is  poorly  defined.  Although  the  data  seem  to  suggest 
that  a nodal  plane  should  lie  south  of  KEV  and  KBS,  and  north  of  NUR  and  GDH,  none  orthogonal 
to  plane  1 could  be  drawn. 

t Pennsylvania  State  University,  University  Park,  PA. 
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TABLE  1-4 

EPICENTRAl  DATA  FOR  THREE  LARGE  EVENTS  OF  1974  SINKIANG  SERIES 


1 

Event 

Dote 

(1974) 

Origin  Time 

Latitude 

(°N) 

Longitude 

CE) 

mb 

Source 

1 

1 1 August 

01:13:55.5 

39.46 

73.83 

D 

EDR 

2 

1 1 August 

20:05:30.9 

39.44 

73.67 

D 

ISC 

3 

1 1 August 

21:21:37.1 

39.46 



73.62 

ISC 

TABLE  1-5 

SOURCE  PARAMETERS  OBTAINED  BY  AMPLITUDE  SPECTRAL  FITTING9 

Event 

Depth 

Dip  Angle 
(deg) 

Slip  Angle 
(deg) 

Strike  of  Nodal  Planes 

1 

10 

90 

0 

75,  165 

2 

10 

90 

0 

30,  120 

3 

6 

50 

-30 

120 

The  discrepancies  between  the  two  types  of  focal-mechanism  solutions  are  larger  than  were 
expected.  Research  into  possible  explanations  for  these  differences,  and  into  the  reliability  of 
the  surface-wave  spectral -amplitude -fitting  procedure.  Is  continuing. 

M.  W.  Shields 


E.  GROUP  VELOCITIES  ACROSS  TIBET  TO  MASHHAD  SRO 

We  are  continuing  our  study  of  group  velocity  dispersion  across  Eurasia,  described  earlier.1 
Although  a large  amount  of  SRO  data  is  now  available,  progress  In  the  inversion  has  been  ham- 
pered by  the  following  factors: 

(1)  The  extreme  tardiness  of  PDE  epicenter  determinations, 

(2)  Closure  of  the  Taiwan  SRO  (TATO)  for  8 months  due  to  technical 
problems, 

(3)  Delays  In  operation  of  Chlengmai  (CHTO), 

(4)  Large  gain  reduction  at  Guam  (GUMO)  which  has  rendered  LP  data 
from  small  events  virtually  unusable. 

We  have,  however,  collected  surface-wave  trains  from  over  200  events  within  and  on  the 
borders  of  Eurasia  and  recorded  at  Mashhad  SRO  (MAIO).  These  include  some  80  events  within 
the  regions  China,  Tibet,  Burma,  and  Mongolia.  The  locations  of  these  events  are  shown 
In  Fig.  1-16.  Many  of  the  paths  to  these  traverse  the  high-altitude  region  consisting  of  the 
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I u .ialuy.is,  1...  i'ibi’Uin  piatoau,  and  ‘.he  Hindu  Kush,  and  the  cross-hutched  area  shown  on 
I ig.  1-16  la  at  altitndes  of  3000  in  or  greater.  Paths  leaving  Mashhad  at  azimuths  of  70‘  to  105 
intersect  this  area  (selected  great-circle  path  shown). 

Figure  1-17  shows  the  envelopes  of  Rayleigh-wave  group  velocities  measured  at  MAIO  for 
events  in  the  azimuth  ranges  9 5"  to  105*  (Southern  Tibet  and  the  Himalayas),  80"  to  90"  (North- 
ern Tibetan  plateau),  and  65"  to  75"  (Hindu  Kush  and  Slnkiang).  respectively.  The  lowest  group 
velocities  are  those  through  Northern  Tibet  (80*  to  90*  azimuth),  and  the  highest  (though  still 
remarkably  low  by  comparison  with  dispersion  for  platform  or  shield  regions)  are  those  which 
cross  the  Himalayas.  Perhaps  the  most  remarkable  feature  of  the  data,  however,  is  the  con- 
tinuation of  these  low  group  velocities  In  the  most  northerly  azimuth  sector,  indicating  that  the 
anomalously  low-velocity  structure  persists  to  the  north  of  the  Tibetan  plateau  as  outlined  by 
the  cross-hatched  region.  These  low  group  velocities  abruptly  cease  for  paths  at  azimuths  of 
less  than  65*.  as  can  be  seen  from  Fig.  1-18  which  shows  dispersion  for  paths  to  Northern  Sln- 
klang  and  Mongolia  compared  with  those  for  Tibet.  (The  shaded  region  here  is  the  envelope  of 
the  group  velocity  regions  shown  in  Fig.  1-17.) 

Perhaps  the  most  striking  illustration  of  the  effect  of  the  Himalayas  and  Tibet  upon  group 
velocity  is  given  by  Fig.  1-19  which  shows  group  velocities  to  events  in  Burma  and  Szechwan. 
Dispersion  to  event  45  (Szechwan,  path  almost  entirely  through  Tibet)  is  up  to  0.5  km/sec  slower 
than  that  to  event  9 (Southern  Burma),  and  it  can  be  seen  that  the  group  velocities  increase 
steadily  as  we  move  the  source  from  north  to  south.  Also  shown  on  the  figure  is  the  dispersion 
calculated  by  Teng15  for  a modal  designed  to  fit  his  data  for  paths  crossing  Southern  Tibet.  The 
paths  he  used  are  similar  to  those  for  event  41  and  the  agreement  is  good.  The  remarkably  slow 
results  for  event  45  are  not  therefore  too  surprising,  since  its  path  comprises  more  than  twice 
as  much  of  the  Tibetan  plateau  as  that  for  event  41. 

Station  CHTO  has  just  come  into  operation  and  as  soon  as  data  are  available  we  shall  be 
able  to  further  determine  the  precise  outline  of  the  low-velocity  region  in  Central  Asia. 

R.G.  North 


F.  GROUP  VELOCITIES  IN  THE  WEST  PACIFIC 

During  the  past  six  n onths,  the  frequency-filter  method  has  been  used  to  obtain  dispersion 
curves  for  over  100  paths  originating  in  the  vicinity  of  the  West  Pacific  and  recorded  at  three 
SRO  stations  - GUMO  (Guam),  SNZO  (New  Zealand),  and  NWAO  (Australia).  Paths  are  shown 
in  Fig.  1-20,  along  with  major  topological  features  of  the  area.  The  region  is  noted  for  pecu- 
liarities in  both  mantle  and  crustal  properties.  In  particular,  the  Fiji  Plateau  is  noted  for  high 
seismic  attenuation  and  slow  P-wave  velocities,14,15  and  the  Ontong-Java  Rise,  to  the  north  of 
the  Solomon  Islands,  has  been  recently  shown  by  refraction  studies  to  have  a crustal  thickness 
appropriate  for  continental  structures.1^ 

Group  velocities  are  shown  for  several  paths  to  GUMO  in  Fig.  I-21(a-e).  Velocities  for 
paths  from  Melanesia  to  Guam  changed  in  a systematic  fashion  as  paths  included  more  of  the 
Fiji  Plateau  and  the  Ontong-Java  Rise.  Paths  from  Samoa  [Fig.  I-21(a)]  gave  the  highest  group 
velocities;  these  are  in  close  agreement  with  velocities  derived  from  the  Dorman  oceanic 
modal  3099  (Ref.  17),  shown  in  the  figure  for  comparison.  Paths  from  Kermadec  [Fig.  I-2i(b)], 
which  crossed  both  the  Fiji  and  the  Ontong-Java  plateaus,  showed  a marked  drop  in  group  veloc- 
ity at  all  frequencies  relative  to  ths  Samoa  paths.  The  results  are  mere  closely  matched  by 
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model  ARC-1  (Ref.  18),  derived  from  a path  from  the  Banda  Sea  to  Japan.  The  lowered  veloc- 
ities and  the  greater  dispersion  at  short  periods  reflect  the  anomalous  mantle  of  the  FIJI  Pla- 
teau and  the  thickened  crust  in  the  Ontong-Java  region.  Paths  from  the  New  Hebrides,  Loyalty 
Islands,  Solomon  Islands,  and  Santa  Cruz  (Figs.  I-21(c),  (d),  and  (e)J  show  even  greater  disper- 
sion at  the  low  frequencies  than  do  the  Kermadec  paths,  reflecting  the  fact  that  a larger  percent- 
age of  the  path  is  through  the  thickened  crustal  region. 


TABL£  1-6 

LAYERED  MODEL  USED  TO  OBTAIN  VELOCITIES 
SHOWN  IN  FIGS.  1-21(0),  (d),  AND  (e) 


Model  A 

Model  B 

h 

a 

p 

P 

h 

a 

P 

P 

2.0 

1,52 

0 

1.03 

25.0 

1.52 

0 

1.03 

1.5 

2.6 

1.09 

2.3 

1.5 

2.60 

0.84 

2.3 

3.0 

5.4 

2.96 

2.6 

3.0 

5.4 

3.0 

2.6 

6.0 

6.1 

3.53 

2.8 

5.0 

6.1 

3.45 

2.8 

17.0 

6.9 

3.98 

2.9 

15.0 

6.9 

3.93 

2.9 

50.0 

8.2 

4.48 

3.3 

50.0 

8.2 

4.46 

3.5 

50.0 

7.9 

4.25 

3.9 

50.0 

7.9 

4.28 

3.4 

250  0 

8.2 

4.45 

3.5 

250.0 

8.2 

4.46 

3.5 

250,0 

9.8 

5.40 

3.96 

250.0 

9.8 

5.40 

3.96 

150.0 

10.0 

5.90 

4.21 

150.0 

10.0 

5.90 

4.21 

- 

11.5 

6.38 

4.95 

- 

11.5 

6.38 

4.95 

The  points  plotted  on  Figs.  I-21(c),  (d),  and  (e)  are  for  velocities  computed  from  two  flat- 
layered models  obtained  using  the  Harkrlder-Rodl  FORTRAN  Inverse  program.  The  models, 
given  In  Table  1-6,  Include  a crustal  thickness  of  27. 5 and  24.5  km,  far  too  high  for  both  oceanic 
and  Island  arc  structures.  Hence,  the  data  support  the  existence  of  an  unusually  thick  crust  in 
the  Ontong-Java  Rise.1*” 

Paths  to  NWAO  from  adjacent  events  in  the  Fiji  Plateau  region  gave  consistent  results,  and 

a systematic  decrease  In  dispersion  with  increasing  oceanic  component  was  observed,  as  shown 

in  Fig.  1-22.  Paths  from  the  Fiji  Plateau  to  SNZO  gave  results  intermediate  between  continental 

and  oceanic,  appropriate  for  island  arc  structures.  Paths  crossing  Indonesia  showed  strong  dls- 

19 

persion.  In  agreement  with  previous  results. 

S.  Seneff 
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Fig.  1-1.  Filter  response  of  elllptlclty  filters  f * 1 - e 
and  g «=  (i  + (e/eo)n)“1//2  described  in  text. 
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Fig.  1-2.  (a)  Rotated  Z and  R LP  component*  from  an  Eastern  Kaaakb  explosion  mw, 
of  S.8  recorded  at  MAIO  SRO.  (b)  Traces  of  (a)  after  polarisation  using  a 3-pole  elllp- 
ticity  filter  g.  (c)  Trace*  of  (a)  after  using  a 6-pole  elliptic ity  filter  g. 
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Fig.  1-3.  (a)  Rotated  Z and  R LP  component*  from  an  Eaetern  Kazakh  exploaion  mb 
of  5.1  obscured  by  an  earthquake  coda,  (b)  Traces  of  (a)  after  polarisation  using  a 
3-pole  elllptictty  filter  g.  (c)  Traces  of  (a)  after  using  a 6-pole  elliptlclty  filter  g. 
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Fig.  1-4.  (a)  Rotated  Z and  R LP  components  from  two  Eastern  Kasakh  explosions  mb 
of  5.1  and  5.3  which  occurred  5 min.  apart.  They  are  obscured  by  an  earthquake  coda, 
(b)  Traces  of  (a)  after  polarisation  using  a 3-pole  elliptic ity  filter  g.  (o)  Traces  of 
(a)  after  using  a 6-pole  elllptlolty  filter  g. 
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Fig.  1-6.  M,  meaaured  from  vertical  Rayleigh  waves  (Mg*)  vs  that 
measured  from  Love  waves  (Mgh. 


Fig.  1-9.  AI.I’A  center  sensor  seismograms  from  five  events  used  in  this  investigation. 
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I'lg.  1-10.  Phase  velocity  dispersion 
data,  estimated  errors,  and  disper- 
sion from  be st -fitting  model. 
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Fig.  1-18.  Group  velocities  to  events  at  azimuths  <65°  (Mongolia. 
N.  Sinkiang).  Also  shown  is  a summary  of  Tibetan  group  veloci> 
ties  (from  Fig. 1-17). 


Fig.  1-19.  Group  velocities  to  events  in  Burma  and  Szechwan, 
together  with  theoretical  dispersion  for  Tibet  given  by  Teng. 
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Fig . 1-20.  Paths  for  which  surface- wave-dispersion  characteristics  have  been  obtained. 


PERIOD  (lid 


PERIOD  UmI 


PERIOD  (Mil 


Fig.  1-21.  Diaper* Ion  curves  for  several  paths  to  Gli MO:  (a)  from  Samoa, 
(b)  from  Kermadec,  (c)  from  Loyalty  Islands,  (d)  from  Fiji  Plateau,  and 
(e)  from  Solomon  Islands. 


11.  SEISMIC  STAKING  AND  BODY-WAVE  MAGNITUDE 


""b 

A.  MOMENTS  AND  STRESS  DROPS  FOR  THE  RAT  1 SI  .AND  SEQUENCE  OF  1965 

tn  a previous  study  of  the  Rat  Island  sequence  of  1965  (Ref.  1),  It  was  shown  that  the  trend  of 

log  (A/T)  or  m.  vs  pulse-period  T,  reported  by  VELA  stations,  favored  an  w-cube  displacement 

0 2 3 4 

model  suggested  by  Savage  rather  than  w-squared  sources  such  as  Brune's  model  * If  geo- 
metrical similarity  was  assumed  for  all  earthquakes.  This  implied  constant  stress  for  all 
earthquakes  in  the  model  with  amplitude  variations  being  controlled  by  the  DC  spectral  level 
which  was  set  equal  to  fj3,  the  inverse  corner  frequency  cubed. 

We  have  now  used  the  empirical  scaling  of  the  synthetic  data  to  estimate  DC  level  of  the 
source  spectra,  and  from  that  the  moment  and  stress  drop  of  the  Rat  Island  data. 

Figure  Il-l(a-b)  shows  the  scatter  of  the  Rat  Island  data  used  to  define  the  trend  of  mb  vs 
period  T.  Of  the  868  aftershocks  selected,  a subset  of  245  events  was  reported  at  all  five  VELA 
stations.  The  mean  and  standard  deviation  of  both  mb  and  T for  each  event  are  shown  in 
Fig.  Il-l(a).  We  see  that  the  trend  of  mb  vs  T is  well  defined  for  the  shorter  periods,  but  seems 
to  deteriorate  for  periods  greater  than  1.2  sec. 

To  quantify  this  trend,  a weighted  linear  regression  line  waB  fit  to  the  data  using  an  itera- 
tive technique  described  by  York.  This  method  allows  individual  weighting  of  every  data  point 
in  both  coordinates.  Following  York's  suggestion,  we  chose  weights  which  are  reciprocals  of 
the  data  variance  in  each  coordinate.  The  resulting  leaM-squares  line  is 

mb  = (3.24  ± 0.08)  + (1.32  ± 0.09)T 
which  is  shown  in  Fig. Ii-l(a). 

Figure  H-l (b)  contains  a histogram  of  average  periods  for  this  data  set  grouped  in  bins 
0.1  sec  wide.  The  distribution  of  periods  is  asymmetric,  peaked  at  0.8  sec  with  90  percent  of 
the  data  falling  In  the  range  G.6  to  1.3  sec. 

Figure  n-2  shows  the  synthetic  curves  of  n^  vs  T for  the  w-cube  model  of  Savage2  as  a 
function  of  comer  frequency  fQ  and  attenuation  parameter  t*.  The  regression  line  of  Fig.  n-l(a) 
is  displayed  as  a pair  of  solid  lines  shifted  by  0.5  m^  unit  to  show  the  effect  of  shifting  the 
theoretical  curves  by  the  same  amount.  The  data  and  synthetic  trends  are  not  perfectly  matched, 
but  are  much  more  similar  than  when  w -square  source  models  are  used.1  The  variation  in  t* 
for  this  empirical  fit  of  model  to  data  is  from  -0.25  < t » < 0.75, 

We  shall  pick  a tabulated  point  on  the  synthetic  mb  curves  in  Fig.  1 1-2  and  first  obtain  the 
DC  spectral  level,  from  which  moment  can  be  estimated.  Under  the  assumptions  of  our  model, 
the  stress  drop  is  constant  for  all  points  along  each  solid  curve  as  fQ  varies. 

The  synthetic  point  in  Fig.  II- 2 with  a corner  frequency  fQ:  I and  t*  = 0.5  has  an  empirical 

of  4.6  on  the  left-hand  scale  and  an  observed  pulse  period  T - 0.90  sec.  This  point  is  a 
little  above  the  lower  Rat  Island  data  trend.  Using  the  Gutenberg  Q factor  for  a distance  cf  60*, 
which  la  about  the  average  VELA  station  distance,  we  obtain  log  (A/T)  “ 1.0  or 

A/T  « 10  mp/aec 

This  value  is  86  times  the  unsealed  synthetic  A/T  value  which  had  a DC  level  set  equal  to  fj3  = 1. 
Thus,  the  DC  spectral  level  of  the  empirically  scaled  displacement  source  in  Fig.  IT-2  is 

_7 

no  = 86  x 10  cm-sec  . 
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Hanks  and  Wyss  proposed  an  adaptation  for  P waves  of  Keilis-Borok's  formula  for  seismic 
moment  which  is 

= no2irpo3  R/r 

where  p = density,  a * P-wave  velocity  in  the  source  region,  R is  the  geometric  divergence 
correction,  and  r is  the  radiation  pattern.  This  includes  a factor  of  0.5  to  account  for  the  free 
surface  reflection  at  the  seismometer.  Hanks  and  Wyss  note  that  SP  amplitudes  are  not  well 
correlated  with  radiation  patterns  based  on  LP  P waves,  and  we  have  no  fault-plane  solutions 
for  there  small  events.  Under  this  circumstance,  we  set  r s 1 although  r is  probably  <1  based 
on  the  fault-plane  solutions  for  large  shocks  of  the  Rat  Island  sequence  obtained  by  Stauder.7 

Taking  as  typical  values  p * 3.0  g/cm3,  a = 6.0  km/sec,  and  R <■  1,9  x 10*  km  at  60* 

(Ref.  8),  we  obtain  for  the  moment 


Mq  “ 6.65  x iO22  dyn-cm 


(11-1) 


Using  this  moment  and  the  corner  frequency  f = 1,  we  can  evaluate  the  stress  drop  for 


this  family  of  curves  by  Brune's  formula,  i.e., 

-2 


7 Mo 
A<r  ■ — y dyn-cm 


where  a is  the  radius  of  the  circular  source  area  given  by 

2.34<; 

*o 


Wf  ta 


Substituting  Eqs.  (11-3)  and  (II-l)  into  (11-2),  we  obtain  the  stress  drop 
■ 25,4  bars  . 


(11-2) 


(II- 3) 


(II-4) 


This  value  of  stress  drop  Inferred  for  the  Rat  Island  earthquakes  with  4 < < 5 agrees  with 

stress  drops  estimated  for  larger  Aleutian  earthquakes  using  spectral  techniques.  Wyss  gt  a^ 
used  LP  P-wave  spectra  to  calculate  stress  drops  of  10  to  47  bars  for  three  Aleutian  earthquakes, 
m^  > 6.0  near  Amchltka  Islands.  The  moment  value  in  Eq.  (II- 1)  is  quite  reasonable  for  - 4.6 
earthquakes  as  extrapolated  down  from  moments  of  larger  earthquakes,8  or  inferred  from 
moment  vs  mL  curves  for  local  earthquakes.6 

Therefore,  it  appears  that  the  use  of  simple  time-domain  models  can  yield  reasonable  mo- 
ments and  stress  drops  for  moderate-slsed  telesetsmtc  earthquakes,  without  resorting  to  LP 
spectral  analysis.  It  is  fortunate  that  DC  levels  Inferred  from  SP  amplitudes  seem  to  agree  with 
actual  LP  measurements.  However,  excellent  bulletin  data  are  required  for  many  events  from  a 
single-source  region,  in  order  to  define  a reliable  trend  of  amplitudes  vs  periods.  Statistically, 
any  single  earthquake  may,  in  tact,  not  fit  the  average  trend  at  all. 

C.W.  Frasier 
R.  Q.  North 


B.  INCORRECT  SEISMIC  SCALING  FROM 

Body-wave  magnitude  (m^)  Is  often  used  to  Indicate  the  spectral  level  of  earthquake  sources 
at  1 Hz  (Refs.  11  and  12).  The  basic  assumption  is  that  an  difference  of  1.0  unit  Implies  a 
factor  of  10  between  the  spectral  levels  of  two  similar  earthquakes  at  a 1-sec  period. 


For  stochastic  source  models,**  this  assumption  is  probably  wrong  because  the  power  level 
at  a single  frequency  cannot  be  calculated  from  a single  amplitude  measurement  of  a complex 
seismogram.  For  deterministic  .source  models,  we  can  clearly  demonstrate  that  is  not  a 
measurement  of  1-Hz  energy  and  that  a significant  error,  which  depends  on  the  earthquake 
source  function  and  attenuation,  results  from  such  an  interpretation. 

To  demonstrate  this,  synthetic  records  for  two  seismic  source  models  were  calculated  as 
recorded  by  four  SP  recording  systems.  Figure  11-3  shows  the  displacement  response,  nor- 
malized to  1 at  i Hz,  of  four  SP  recording  systems  - WWSSN,  VELA,  SRO,  and  LASA.  The 
SRO  response  was  calculated  from  the  system  response  described  in  Sec.  IV-A  of  this  SATS. 

Both  SRO  and  LASA  systems  are  digitally  recorded  with  anti-alias  filters  and  no  galvonom- 
eter,  whereas  the  VELA  and  WWSSN  systems  have  coupled  seismometer  and  galvonometer 
combinations.**'*4 

Figure  II-4  displays  the  impulse  displacement  response  of  the  four  systems,  including  a 
small  amount  of  causal  earth  attenuation  t*  T/Q  = 0,1,  where  T equals  the  travel  time  and 
Q the  average  quality  factor  for  the  teleseismlc  path.  The  ripple  in  the  coda  of  LASA  and  SRO 
responses  is  due  to  the  sharp  cutoff  of  the  anti-alias  filter  in  each  system,  the  SRO  filter  being 
especially  sharp.  One  also  sees  a large  group  delay  for  the  SRO  response  (see  Sec.  IV-A). 

The  impulse  responses  in  Fig.  11-4  were  further  attenuated  and  convolved  with  an  w-aquare 
source  by  Brune*  and  an  to -cube  source  by  Savage*  to  produce  a suite  of  synthetic  records  for 
each  recording  system.  Figure  II- 5 displays  synthetic  records  for  the  co-square  source  with  a 
comer  frequency  of  f * 0.5  Hz  and  t*  = 0.5.  Although  the  system  characteristics  in  Flgs.II-3 
and  11-4  vary  considerably,  the  overwhelming  effect  of  earth  attenuation  and  source  decay  above 
the  corner  frequency  produces  similar-looking  records  for  all  four  seismometer  systems. 

Synthetic  records  were  calculated  for  a range  of  corner  frequencies  f.  from  2*  = 8 Hz  to 
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2 a 0.00391  Hz  for  each  source  model  assuming  two  t*  values,  0.5  and  1.0.  Relative  magni- 
tudes were  computed  for  each  seismogram  by  calculating  log(A/T)  where  A is  peak  amplitude 
and  T la  peak  period,  and  A is  corrected  for  seismometer  response  at  period  T.  Due  to  ob- 
served dispersion  In  the  waveforms  and  ambiguities  in  the  choice  of  peaks  used  for  A and  T by 
an  analyst,  A was  measured  as  the  peak  value  of  the  envelope  of  the  trace,  and  T was  calculated 
from  the  instantaneous  frequency  f = l/T  at  that  envelope  peak.  This  method  produces  magni- 
tudes very  close  to  the  values  obtained  by  an  analyst,  and  produces  more  continuous  T 
measurements. 

In  Fig.  II-6(a),  the  difference  between  log(A/T)  and  log  of  the  source  at  1 Hz  is  displayed 
for  the  Brune  source  as  a function  of  comer  frequency  fQ.  The  log  (A/T)  values  are  averages 
of  all  four  recording  systems.  The  standard  deviations  of  the  log  (A/T)  distribution  is  very 
small  (<0.03)  for  all  corner  frequencies,  so  that  relative  magnitudes  are  practically  identical 
for  all  four  recording  systems  in  spite  of  the  variations  in  displacement  response.  Figure  II-6(b) 
shows  similar  plots  for  the  w-cube  Savage  source. 

Although  the  extremes  of  corner  frequencies  are  outside  the  range  expected  for  teleseismlc 
data,  the  region  of  most  rapid  change  in  the  curves  is  between  2.0  and  0.0625  Hz,  which  includes 
most  teleseismically  recorded  magnitudes.  Bars  in  Figs.ll-6(a)  and  (b)  indicate  the  range  of 
comer  frequencies  expected  for  m^  4.5  and  6.0  as  the  stress  drop  varies  from  i to  100  bars 
for  each  magnitude,  calculated  by  North  and  Frasier  (see  Sec.C  below).  Thus,  he  rapid  change 
in  bias  of  log  (A/T)  about  log  of  the  source  at  1 Hz  cannot  be  adequately  removed  oy  a constant 
shift  of  log  (A/T)  such  as  the  Gutenberg  Q factor. 


The  basic  cause  of  this  non-constant  bias  is  the  variation  in  observed  periods  as  a function 
of  corner  frequency.  Figures  II-7(a)  and  (b)  show  the  average  pulse  period  for  each  source  for 
the  two  t*  values.  The  periods  are  averaged  over  the  four  recording  systems,  and  the  period 
distributions  have  standard  deviations  e.f  <0.08  sec  for  the  Prune  source  at  all  corner  frequen- 
cies. For  the  Savage  source,  the  period  distributions  have  standard  deviations  <0.11  sec  for 
t*  - 0.5,  and  <0.20  sec  for  t*  = 1.0,  the  largest  deviations  occurring  at  corner  frequencies 
<0.0625, 

Using  the  observed  period,  we  can  calculate  the  expected  frequency- domain  amplitude  at 
that  period  by  taking  the  product  of  source  function  and  attenuation  and  dividing  by  T.  A com- 
parison of  this  frequency-domain  estimate  of  log  (A/T)  with  the  observed  log(A/T)  on  the  syn- 
thetics is  shown  in  Figs.  11- 8(a)  and  (b).  The  agreement  is  nearly  perfect  and  indicates  that 
observed  time-domain  magnitudes  for  deterministic  sources  predict  almost  exactly  the  expected 
spectral  amplitudes  at  the  pulse  period  If  the  source  function  and  attenuation  are  included. 

From  this  synthetic  record  experiment,  it  is  clear  that  m^  variations  cannot  be  used  to  infer 
seismic  source  scaling  at  1 Hz  without  detailed  knowledge  of  the  source  model.  Including  corner 
frequency  and  earth  attenuation  along  the  teleseismic  path. 

C.  W.  Frasier 

C.  mb-FREQUENCY  RELATIONS  FROM  SEISMIC  SCALING  LAWS 

In  an  earlier  attempt1  to  reproduce  a linear  trend  of  log  (A/T)  vb  T (where  A is  the  maxi- 
mum peak-to-peak  amplitude,  and  T the  dominant  period)  of  P waves  observed  for  aftershocks 
of  the  Rat  Island's  earthquake  of  1965  we  have  generated  synthetic  seismograms  for  a range  of 
corner  frequenciea  tQ  for  the  Brune  (w2  far- field  spectral  decay  below  the  corner  frequency) 
and  Savage  (w3)  source  models.  We  found  that  the  data  trend  of  log  (A/T)  vs  T was  best  fitted 
by  the  w3  model  with  t*  = 0.5.  An  inherent  assumption  in  the  generation  of  the  synthetic  seis- 
mograms was  that  of  geometrical  similarity.  Implying  constant  stress  drop,  l.e., 

logMQ  » -3  logfQ  + C <M0  * seismic  moment)  (11-5) 

where  the  constant  C depends  upon  the  stress  drop  Ao. 
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Figure  11-9  ehows  the  variation  of  with  log  (fQ)  for  the  w and  w source  models,  for 
t*  = 0.5  and  1.0,  For  a given  corner  frequency  f0,  the  difference  in  mfe  given  by  the  two  values 
of  t*  la  ~0,7  unit  at  fQ  = 8 Hz,  decreasing  to  <0.2  unit  at  the  highest  values  of  f0<  This  variation 
is  caused  by  the  increase  in  dominant  period  of  the  P wave  as  f decreases.  It  can  be  seen,  as 
remarked  previously,1  that  the  u>2  model  predicts  higher  values  of  m^  for  a given  f than  the 
w3  model  for  f < 0.5  Hz, 

We  show  elsewhere  (see  Sec.  B above)  that  the  stress  drop  Implied  by  the  fit  of  the  Rat  Is- 
land data  to  the  w2  model  la  25  bars,  derived  by  correcting  the  spectrum  for  a source  with 
fQ  = 1 Hz  back  to  the  source  to  obtain  a moment  of  6.7  x iO22  dyn-cm.  Using  this  result,  we 
may  compute  the  precise  form  of  C (Ao)  In  Eq.  (II-5). 

By  combining  Fig.  11-9  and  relation  (II-5),  we  may  obtain  the  variation  of  log(MQ)  with  mb 
for  the  u2  and  w3  models  for  any  particular  values  of  Act  and  t* . Figure  tt-10  shows  moment 
as  a function  of  m^  for  Ao  » 10  bars  and  t*  * 0.5.  Different  values  of  stress  drop  A<r  merely 
shift  the  curves  up  or  down;  changes  In  t*  cause  variations  in  m^  for  a given  UQ  which  are 
small  for  tn^  > 5.0, 
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1i  can  ho  seen  from  Fig.  11-10  that  the  most  striking  difference  between  the  a;2  and  w3  model 

is  that  the  latter  predicts  an  tipper  limit  to  ~ 6.0,  whereas  the  former  predicts  that  will 

increase  with  moment  without  limit.  For  the  smallest  stress  drop  we  have  considered  (1  bar), 
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the  u-  model  predicts  an  m.  - 7 . f.  for  M 10  dyn-cm.  Such  large  values  of  m,  have  never 

n o D 

been  observed,  even  though  events  of  this  moment  appear  to  occur  about  once  every  i to  2 years. 

A considerably  more-sophisticated  source  model  than  the  simple  ones  we  have  considered, 
incorporating  not  only  the  effects  of  source  dimension  but  also  fault  width  and  rise  time,  has 
been  described  by  Geller.17  The  high-frequency  asymptote  of  far-field  spectrum  produced  by 
tids  model  is  ic"1  decay,  but  it  also  Includes  regions  of  w~2  and  w*1  decay.  From  the  relations 
given  by  Geller  between  Mq  and  Mg,  and  Mg  and  m^,  we  may  deduce  the  variation  of  Mq  with 
m^,  and  this  also  is  shown  in  Fig.  11-40.  Geller  assumed  that  Is  dtrectly  related  to  the 
spectral  density  at  1 Hz;  this  is  not  strictly  correct,  as  he  mentions  and  ao  is  also  shown  else- 
where In  this  report  (Sec.R  above).  However,  it  can  be  seen  that  the  variation  of  log(M„)  with 
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mfe  that  he  obtains  is  very  similar  to  that  produced  by  the  much  simpler  w model  we  have  con- 
sidered, and  also  implies  an  upper  limit  to  m^  of  ~6.0. 

We  have  no  data  on  the  variation  of  MQ  with  against  which  to  test  the  model  predictions, 
but  there  does  exist  a large  amount  of  data  on  the  rate  of  occurrence  of  events  of  a particular 
mb  (known  as  the  m^- frequency  curve).  It  has  also  been  shown  that  for  logMQ  > 25,  the 
frequency  curve  is  remarkably  linear  and  given  by 

log  N n 17.47  - 0.61  log  Mq  . (II-6) 

Evidence  that  this  linear  relation  persists  for  smaller  events  is  provided  by  data  from  Southern 
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California  where 

log  M = i,5  M.  +16.0  , for  20  < log  M < 25 

O L O 

(M.  - local  magnitude)  and  the  M.  -frequency  curve  is  remarkably  linear  for  1 < M.  <7  with  a 
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slope  very  close  to  unity.  Thus,  for  Southern  California,  at  least, 

logN  - constant  - 0.67  logMo  . (II-7) 


The  slopes  of  Eqs.  (II- 6 ) and  (II-7)  are  very  close,  and  we  therefore  assume  that  Eq.  (II-7)  holds 
for  20  < log(M0)  < 31. 

By  combining  Fig.  11-10  with  relation  (11-6),  we  may  obtain  m^-frequency  curves  for  the 
source  models  we  have  considered.  Those  for  Act  - 50  bars  and  t*  =0.5  are  shown  in  Fig.  11-11. 
Changes  in  Act  will  shift  the  curves  up  or  down;  changes  in  t*  produce  appreciable  differences 

only  for  m.  < 5.0.  Also  shown  on  this  figure  is  an  m. -frequency  curve  derived  from  the  ISC 

20  D D 
dataset  for  1964  to  1973,  Below  m^  - 5,  the  shape  of  this  curve  is  primarily  dictated  by 

network  detection  characteristics  and  bearB  no  known  relation  to  the  actual  occurrence  of  events 

3 

with  a particular  m^.  It  can  be  seen  that  the  data  lie  between  the  curves  produced  by  the  u and 

w2  models.  No  possible  combinations  of  stress  drop  and/or  t*  can  reproduce  the  ISC  data 

curve.  Our  theoretical  curves  are,  of  course,  critically  dependent  upon  the  slope  of  Eq.  (II-6), 

but  unreasonably  large  changes  in  the  slope  are  required  to  make  either  model  match  the  data. 

5 /2 

It  is  interesting  to  note  that  an  w ' scaling  law  would  produce  an  m. -frequency  curve  lying  be- 

23  ° 24  22 

tween  those  for  u and  w ; such  a scaling  law  has,  in  fact,  been  suggested  in  the  literature.  ' 

It  can  be  seen  that,  provided  a linear  relationship  of  the  form  Eq.  (11-6)  holds,  none  of  the 

scaling  laws  mentioned  here  (w2,  w3,  w5/2,  or  Gellers  model)  predicts  that  the  m^-frequency 
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<•  lit- vt*  is  linear  except  at  vi  ry  small  magnitudes.  Thus,  "b-valuea"  are  essentially  meaningless 

2 } 

except  over  very  limited  magnitude  ranges.  M.  Wyss  has  claimed  that  low  b-values  Indicate 
high  stress  in  the  source  region  and  vice  versa,  assuming  a relation  ot  the  form  of  Eq.(II-6). 
Our  results  show  that  the  only  effect  of  changing  3tress  drop  is  to  move  the  m^-frequency  curve 
up  or  down;  for  any  one  of  the  source  models  postulated,  the  shape  of  the  curve  is  completely 
independent  of  stress  drop.  R North 

C.  W.  Frasier 


O.  SATURATION  OF  THE  mfa  SCALE 

The  study  described  above  demonstrates  that  information  about  the  seismic  source  can  be 
obtained  from  data  on  the  saturation  of  the  scale  at  large  magnitudes.  The  only  empirical 
data  of  this  kind  that  are  available  are  contained  in  global  earthquake  catalogs  such  as  the  Bul- 
letin of  the  International  Seismological  Center  (ISC),  and  it  is  well  known  that  there  are  problems 
in  assigning  rnb  values  to  large  events. 

The  investigation  of  saturation  of  the  m^  scale  Involves  estimating  the  shape  of  the  frequency 
magnitude  relationship  for  large  events.  Previous  studies*  indicated  that  network  estimates  of 
mb  for  large  events  will  be  affected  by  the  detection  probability  curves  for  individual  stations 
in  the  network,  which  will  fall  off  due  to  instrument  saturation  or  clipping.  Tne  result  of  clipping 
will  be  to  Introduce  a negative  bias  into  network  ni^  estimates  for  large  events,  in  exactly  the 
same  way  that  detection  thresholds  introduce  positive  biases  at  low  m^  values.  The  effect  of 
clipping  is  therefore  in  the  same  direction  as  saturation  due  to  the  seismic  source,  and  it  is 
not  easy  to  separate  these  two  effects. 

If  the  mb  values  listed  in  the  ISC  Bulletin  are  taken  at  face  value,  saturation  appears  to 
occur  in  the  vicinity  of  mfa  = 6.5.  Figure  II- 1 2(a)  shows  the  frequency-m^  values  for  1966-70, 
for  all  regions  of  the  world.  Notice  the  large  apparent  "b-value,"  which  is  somewhat  subjective 
but  which  might  be  as  large  as  1.5.  Selection  of  a smaller,  more  well-controlled  network  of  28 

stations  with  good  reporting  characteristics  leads  to  essentially  the  same  result  [Fig.  1I-12(b)|. 
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Inclusion  of  station  biases,  as  determined  by  North,  does  not  change  the  curve  significantly. 

A very  substantial  change  is  observed,  however,  if  the  VELA  arrays  located  in  the  Western 
U.  S.  are  used  alone.  Figure  11-13  shows  frequency-m^  plots  for  the  four  arrays  BMO,  TFO, 
UBO,  and  WMO,  with  the  station  biases  Included.  The  agreement  between  stations  is  excellent, 
and  the  mean  frequency-m^  curve  appears  to  have  a shallow  slope  (0.93)  at  low  magnitudes  which 
begins  to  depart  from  linearity  at  about  mb  = 5.7,  and  saturation  Is  indicated  at  > 7.0. 

Clearly,  either  the  m^  values  for  large  eventB  reported  by  the  VELA  arrays  are  overesti- 
mated, or  the  network  average  values  listed  by  the  ISC  are  subject  to  a severe  negative  bias 
for  large  events.  All  evidence  obtained  to  date  suggests  that  the  latter  is  the  case.  The  VELA 
arrays  were  operated  by  trained  personnel,  and  the  instruments  were  carefully  calibrated.  More 
importantly,  they  were  equipped  with  a low  gain  channel  which  gave  them  a wide  dynamic  range. 
Also,  the  form  of  the  seismicity  curve  in  Fig.  11-13  is  what  one  might  expect,  since  it  is  almost 
Identical  to  the  shape  of  the  frequency-Ms  curve  discussed  by  Chtnnory  and  North.16 

Following  the  arguments  given  In  Ref.  16,  and  assuming  that  the  VELA  seismicity  curve  is 
correct,  we  can  make  a first  empirical  estimate  of  the  shape  of  the  mb*  moment  relationship. 

This  is  shown  in  Fig.  11-14.  The  data  are  consistent  with  a linear  relation  between  Mfl  and  m^ 


fur  loss  thnn  5.7  or  S.H.  The  constant  in  this  relationship  Is  not  well  known,  and  we  have 
used  0.5  in  constructing  l‘m.  11-14.  Changing  the  constant  will  move  the  m^- moment 

curve  parallel  to  the  log  (moment)  axis  by  a small  amount. 

The  m. -moment  relation  shown  lies  between  the  a ^ and  models  discussed  in  the  pre- 
ceding contribution  of  this  report,  but  there  is  n suggestion  that  it  is  somewhat  closer  to  the 

M.  A.  Chlnnery 
.1.  C.  Johnston 
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Fig.  II-l.  (a)  Average  magnitude  va  average  period  T for  245  events 
commonly  recorded  at  five  VELA  etatlona.  Error  bare  Indicate  standard 
deviations  of  mjj  and  T populations  about  average  values.  Least* squares 
line  through  mean  values  is  shown,  (b)  Histogram  of  average  periods  in 
(a)  grouped  in  bins  0.1  sec  wide. 


Fig.tt-2.  Family  of  synthetic  record  ampli- 
tudes predicted  by  Savage  source  model, 
shifted  to  fit  trend  of  Rat  Island  data.  Right- 
hand  mi,  scale  refers  tc  upper  heavy  Rat  Is- 
land trend  line;  left  scale  refers  to  lower 
Rat  Island  trend  line.  Hatchured  line  Is 

vs  T for  U.  S.  explosions  (Marshall10),  left- 
hand  m jj  scale.  Synthetic  amplitudes  are 
calculated  for  corner  frequencies  0.0629  to 
8 Hz,  and  t*  * 0.1  to  1.0. 


Fig.  11-3.  Displacement  response  of  four 
seismometer  systems  - WWSSN,  VELA, 
SRO,  and  LASA. 


F»g.  II-6.  Plot*  of  average  [log  (A/T)  - 
log  source  at  1 sec)  for  two  source 
models  and  two  attenuation  parameters. 
Log  (A/T)  values  are  averaged  for  four 
recording  systems,  (a)  Grune  source 
(w-2).  (b)  Savage  source  (w3). 
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Fig.  11-7.  Average  pulse  period  T 
on  synthetic  records  for  two  source 
models  and  two  attenuation  param- 
eters. T values  are  averaged  for 
four  recording  systems,  (a)  Brute 
source  (w*2)j  (b)  Savage  source 


coasts  nmutNCv  fiat 

Fig.  11-8.  Plots  of  average  (log  (A/T)  - 
(log  source  at  T sec  + log  atten  at 
T - logT)]  for  two  source  models  and 
two  attenuation  parameters,  (a)  Brute 
source  («•*)(  (d)  Savage  souxoe  (w-3). 


Fig.  II- 10.  Seismic  moment  me  a func- 
tion of  m|)  for  eource  models  of  Brune 
(w2).  Savage  (w  3),  and  Qeller,  for  Ac  ■ 
10  bare  ana  t*  = 0.5. 


Ftg.II-11.  m^- frequency  curvet  for  u2 

and  u*  models  with  Ac  ■ 50  bare  and 
t*  • 0,5,  Alao  shown  Is  curve  for  ISC 
dataset,  1964  to  19?5.  All  curves  re- 
duced to  annual  occurrence  rates, 
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rig.  11-12.  (•)  Frequency-in^  data  for  all  events  in  ISC  Bulletin  for  1966  to  1970. 

(b)  values  recomputed  using  a 28-station  network,  and  requiring  detection  by 

at  least  three  stations.  Open  circles  show  changes  that  result  when  station  biases 
are  included. 


Fig.  11-13.  Seismicity  for  1966  to  1970  observed  by  four  VELA  arrays 
in  Western  U.S.  Station  biases20  have  been  Included. 


Fig.  11-14.  mv-moment  relationship  de- 
duced from  seismicity  curve  in  Fig.  11-13, 
assuming  moment-frequency  curve  quoted 

by  Chtnnery  and  North16  is  correct.  Em- 
pirical M -moment  relationship  from  same 
16  " 

study  is  shown  for  comparison. 


III.  MISCELLANEOUS  STUDIES 


A.  AN  ANALYSIS  OF  THE  MAIO  RECORDING  OF  THE  20  MARCH  1976 

EASTERN  KAZAKH  EVENT 

The  multiple  event  discussed  in  the  previous  SATS1  was  analyzed  by  deconvolution  and 
maximum  entropy  cepstral  analysis.  We  used  the  SP  SRO  recording  made  at  MAIO,  Although 
SRO  data  are  the  very  best  available  for  digital  analysis  in  terms  of  their  dynamic  range  and 
SNR,  the  choice  of  the  MAIO  recording  led  to  two  immediate  problems.  The  first  concerns  the 
SRO  SP  recording  system  response  which  is  discussed  in  Sec.  IV  of  this  report.  Due  to  the 
extreme  anti-alias  filter,  the  SRO  SP  response  is  characterized  by  substantial  group  delay  and 
S-Hz  ringing.  Clearly,  any  attempt  to  understand  SRO  SP  recordings  should  take  this  effect 
into  account.  The  second  difficulty  1s  due  to  propagation.  Eastern  Kazakh,  being  about  20*  away 
from  MAIO,  radiates  energy  through  one  of  the  principle  discontinuities  In  the  upper  mantle.  The 
result  is  artificial  echoes  in  the  wavetrain. 

Our  approach  to  the  first  difficulty,  that  due  to  the  SRO  SP  recording  system  response,  was 
to  design  an  optimum  least-squares  inverse  filter  which  deconvolves  the  recording  system  lm- 
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pulse  response  out  of  our  seismogram.  This  filter  was  computed  by  a program  due  to  Robinson 
which  also  optimizes  the  time  delay  of  the  resulting  filter  output.  Our  filter  was  41  points  (2  sec) 
long  and  was  computed  in  the  presence  of  5-percent  additive  white  noise.  The  optimum  delay 
time  was  1.4  sec.  Figure  III-l  shows  the  result  of  convolving  this  filter  with  the  SRO  SP  Im- 
pulse response.  As  can  be  seen,  the  Inverse  filter  "spikes"  the  Input  Impulse  response  sharply 
at  1.4  sec  delay.  The  other  smaller  disturbance  at  approximately  4 sec  delay  Is  due  to  the  filter 
seeing  the  offset  at  the  end  of  the  impulse  response. 

The  results  of  applying  this  inverse  filter  to  the  MAIO  recording  of  the  20  March  1976  event 
art  .hown  as  the  first  four  traces  in  Fig.  111-2.  Trace  1 is  the  original  recording.  Trace  2 Is 
the  original  recording  filtered  with  a 3-pole  Butterworth  bandpass  filter  with  low  and  hlgh-3-dB 
points  at  0.5  and  4.0  Hz,  respectively.  Trace  3 is  the  original  recording  deconvolved  with  the 
inverse  filter  described  above.  The  fourth  trace  is  trace  3 filtered  with  the  same  bandpass  filter 
used  to  produce  trace  2. 

Looking  at  trace  3,  one  can  easily  see  that  the  Inverse  filter  removes  the  group  delay  In  the 
SRO  SP  response.  The  filter  delay  of  1.4  sec  has  been  removed  by  shifting  the  trace.  The  effect 
of  the  filter  on  noise  outside  the  signal  band  can  also  be  seen.  This  noise  Is  removed  In  trace  4 
by  bandpass  filtering)  however,  the  phase  of  the  bandpass  filter  Is  Introduced.  So,  trace  4 can 
be  thought  of  as  the  original  recording  equalized  In  amplitude  from  0.5  to  4.0  Hz  with  a 3-pole 
Butterworth  phase  response. 

In  the  first  20  sec  of  the  deconvolved  seismogram  (trace  3),  four  primary  phases,  two  depth 
phases,  and  three  echoes  can  be  Identified.  The  first  arrival  and  Its  associated  depth  phase  and 
echo  are  very  faint.  Its  pP  delay  is  somewhere  in  the  3 to  5 sec  range.  The  second  arrival 
exhibits  a very  clear  pP  phase  with  a delay  of  2.4  sec,  and  an  echo  with  a delay  of  8.8  soc  and 
relative  amplitude  of  0.833.  The  third  arrival  has  no  apparent  pP  phase,  and  its  associated 
echo  Is  mixed  in  with  the  fourth  arrival.  This  last  arrival  has  a clear  echo,  but  Its  pP  phase 
Is  lost  In  the  increasing  coda. 

The  results  of  anti-echo  filtering  the  deconvolved  seismogram  are  shown  as  trace  5.  The 
procedure  consisted  of  computing  the  inverse  e-transform  of  the  echo)  i.e.: 


, b,-i  , b 2 2b  3 3b  . 
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and  convolving  with  the  echoed  seismogram.  The  echoes  EJ#  E2,  and  E4  are  removed  from 
the  trace.  This,  then,  is  the  final  "cleaned-up"  seismogram.  Here,  the  first  motions  of  P2, 
Pj,  and  P4  can  easily  be  seen  to  be  upward.  The  first  motion  of  Pj  cannot  be  identified. 

A maximum  entropy  eepstrum  of  the  deconvolved  seismogram,  trace  4 in  Fig.  1II-2,  is 
shown  in  Fig.  III-3.  The  three  most  prominent  peaks  occur  at  2.6,  4.7,  and  8.9  sec,  respec- 
tively. The  last  can  be  identified  as  the  travel-time  curve  echo  delay.  The  first  two  are  pos- 
sibly pP  times  of  the  primary  arrivals.  In  particular,  4.7  sec  may  correspond  to  the  pP  delay 
of  P4  which  cannot  be  identified  in  the  seismogram. 

On  the  basis  of  this  analysis  we  conclude  that,  in  the  first  20  sec  of  the  MAIO  recording, 
only  one  arrival  appears  to  be  definitely  lacking  a pP  phase.  That  arrival  was  P3  which  ar- 
rived at  approximately  04:08*08.5  on  the  original  seismogram. 

D.  W.  McCowan 
T.  E.  Landers 


B.  CRUSTAL-PHASE  TRAVEL  TIMES  FROM  NTS  EXPLOSIONS: 

PRELIMINARY  RESULTS 

We  have  compiled  an  extensive  data  set  of  travel  times  from  NTS  explosions  recorded  mainly 
at  Long  Range  Seismic  Measurements  (LRSM)  stations  in  the  Western  U.S.  It  was  our  original 
intention  to  rely  on  the  Oeotech  shot  reports  for  tabulated  observations.  However,  we  found  that 
there  were  no  S-phase  times  and  very  few  P*  times  in  these  reports.  So,  in  order  to  acquire 
more  data,  we  reread  the  original  35-mm  LRSM  film  recordings.  A side  benefit  of  this  pro- 
cedure was  observational  uniformity,  since  all  the  data  were  read  by  one  analyst.  We  arbi- 
trarily limited  ourselves  to  explosions  for  whloh  there  were  announced  yields,  giving  us  data 
from  25  events.  A map  of  the  stations  is  Included  as  Fig.  111-4.  Unfortunately,  because  the 
LRSM  sites  continually  changed  location  during  the  recording  period,  there  were  few  common 
stations. 

These  travel-time  observations  were  checked  for  errors  and  erroneous  readings.  They 
were  then  placed  in  our  Data  Analysis  and  Display  System  (DADS)  for  analysis  and  display. 
Altogether,  we  had  171,  76,  54,  47,  183,  and  85  Pg,  Sg,  P*,  S*,  Pn,  and  Sn  observations, 
respectively.  As  can  be  seen,  there  were  fewer  S-phase  observations  than  corresponding 
P-phase  observations  and,  in  both  cases,  more  n-phase  observations  than  either  g or  * ob- 
servations. This  is  not  unusual  and,  most  likely,  is  due  to  the  variability  of  the  crust  over  the 
station  coverage. 

Plots  of  residuals  tor  all  six  phases  as  compared  with  those  from  a Jeffreys -Bullen  (J-B) 
crustal  model  are  shown  in  Figs.  I1I-5  to  III-7.  The  model  is  a two-layer  crust  over  a mantle 
halfspace  with  values  chosen  to  fix  the  J-B  crustal-phase  tables.  Layer  parameters  are  given 
in  Table  III-l.  Least-squares  lines  fitted  to  the  corresponding  residuals  are  also  shown  in 
Figs.m-5  to  IU-7. 

Two  sources  of  error  can  be  easily  seen  in  these  figures.  The  first  is  a random  observa- 
tional error  which  is  Independent  of  distance.  Comparing  Figs.  IQ-5(a),  III-6(a),  and  III-7 (ax) 
with  Figs.  Ill- 5(b),  IH-6(b),  and  III-7(b),  one  can  see  that  the  random  error  component  is  at 
least  twice  as  great  for  S phases  as  for  corresponding  P phases.  We  feel  that  this  error  is 
due  to  inaccurate  time  picks.  The  S phases,  coming  later  in  the  P coda,  are  harder  to  pick 
and  are  therefore  less-accurately  read. 
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TABLE  lll-l 

CRUSTAL  MODEL  FOR  THE  WESTERN 

UNITED  STATES 

Thickness 

j J-B  P 

J-B  S 

P 

S 

Layer 

(km) 

| (km/sec) 

(km/sec) 

(km/sec) 

(km/»ec) 

1 

15 

! 

3.37 

6.03 

3.50 

2 

18 

1 6.50 

3.74 

6.73 

3.80 

3 

00 

! 7.85 

1 

j 

4.42 

7.88 

3.<»7 

The  second  error  component  is  a trend  which  clearly  depends  on  distance.  It  can  be  seen 
best  in  Fig.  III-7(a)  for  the  Pn  phase.  After  removing  the  random  error  due  to  observation,  we 
estimate  that  the  absolute  value  of  the  remaining  error  in  the  Pn  residual  is  approximately 
0.4  sec/deg.  We  feel  that  this  error  component  is  due  to  the  effects  of  crustal  heterogeneity 
in  the  Western  U.S.  For  example,  stations  near  the  NTS  will  only  see  the  Basin  and  Range 
crustal  structure  because  that  is  where  the  NTS  is  located.  Stations  further  away  will  see  other 
crustal  structures,  depending  on  how  fur  away  they  arc  from  the  NTS  and  along  which  azimuth. 

It  is  one  of  our  objectives  to  sort  the  travel  times  by  azimuths  and  attempt  to  Identify  the  effects 
of  varying  crustal  structures. 

The  data  shown  in  Figs.  III-5  to  111-7  can  be  used  to  perform  a refraction  survey  of  the 
Western  U.S.  The  procedure  is  the  same  as  that  commonly  used  in  petroleum  exploration  (for 
example,  see  Dobrin  ).  The  P-  and  S-wave-velocity  results  are  given  in  Table  III-l.  Our 
model  favors  higher  P and  S velocities  in  both  crustal  layers,  a Pn  velocity  very  close  to  the 
J-B  value,  and  a lower  Sn  velocity  than  predicted  by  the  J-B  tables.  The  most  accurate  veloci- 
ties are  Pg  and  Pn.  P*  appears  to  be  variable  between  sites,  and  all  the  S phases  show  con- 
siderable scatter. 

Our  data  were  not  accurate  enough,  before  sorting  into  azimuths,  to  resolve  the  layer 
thicknesses.  The  P*  and  S*  times  Indicate  a first-layer  thickness  of  27  and  14  km,  respectively. 
We  do  not  consider  either  value  significant.  On  the  other  hand,  the  very  small  Intercept  on  the 
Pn  residual  plot  |Flg.  Ill-7(a)l  implied  a Moho  depth  close  to  the  J-B  value  of  33  km.  Since  these 
are  the  data  with  the  least  scatter,  we  do  consider  this  value' significant.  To  summarize,  our 
preliminary  analysis  of  this  crustal-phase  data  indicates  that  the  J-B  crustal  thickness  and  Pn 
velocity  are  good  approximations  to  the  Western  U.S.  The  J-B  crustal  layers  have  lower  P- 
and  S-wave  velocities  than  our  data  suggest. 

These  preliminary  results  generally  agree  with  those  from  other  investigations.  Our  crustal 
thickness,  33  km,  agrees  with  Glover  and  Alexander's4  crustal  thickness  for  the  NTS  as  well  as 
with  Warren  and  Healy's5  contour  map  of  Moho  depth  in  the  Western  U.S.  Our  Pn  velocity, 

7.88  km/sec,  is  the  same  as  that  determined  by  long-profile  refraction  surveying  reported  by 
Eaton.6  In  short,  these  data  support  the  substantial  tectonic  research  done  on  the  Basin  and 
Range  province. 

We  intend  to  use  these  data  in  an  experiment  to  determine  how  much  crustal-phase  travel- 
time  observations  Improve  teleselsmic  hypocenter  locations.  After  fixing  a structure  for  the 
Western  U.S.  on  the  basis  of  those  data,  we  will  start  relocating  NTS  explosions  and  earthquakes 


in  the  general  source  region.  We  expect  that  these  results  will  aid  in  evaluating  the  location 
capabilities  of  tin;  near-in  networks. 

D.  W.  McCowan 

II.  K.  Needham 

C.  THE  DEPTH  OF  SEISMICITY  IN  AN  IRANIAN  AFTERSHOCK  SEQUENCE 

The  lack  of  nearby  stations  and  poorly  known  velocity  structure  of  parts  of  Asia  causes  the 
teleseismically  found  focal  depths  of  earthquakes  in  these  regions  to  be  badly  controlled  — an 
effect  seen  clearly  in  the  Zagros  mountains  of  Iran,  where  the  reported  deep  events  are  nearly 
all  poorly  recorded  (Fig.  IH-8). 

The  Gir  (1972)  aftershock  sequence  in  Southern  Iran  was  relocated  by  the  master-event 
technique7  using  ISC  Bulletin  data,  to  try  and  precisely  determine  depths  relative  to  the  master; 
the  advantage  of  this  method  is  that  the  dependence  of  computed  depths  on  the  assumed  earth 
model  is  reduced. 

The  main  destructive  event  (magnitude  6.0)  was  a thrust  preceded  by  a small  foreshock, 
making  onsets  difficult  to  read;  but  in  the  absence  of  a large  event  as  an  alternative,  it  was 
retained  as  the  master.  The  relocations  are  shown  in  Figs.  III-9(a)  and  (b)  and  have  the  same 
strike  as  the  fault-plane  solution  for  the  master;  they  also  define  a plane  dipping  NE  at  60*  to 
70*  - an  effect  dramatically  seen  when  their  position  vectors  are  plotted  on  a stereographic 
projection  centered  on  the  master  (Fig.  IH-9(c)].  Locations  are  all  relative  to  the  master  event, 
whose  initially  chosen  position  determines  the  absolute  locations.  Great  care  was  taken  to  ex- 
clude bad  data  from  the  relocations;  all  records  from  WWSSN  stations  closer  than  30*  were 
reread  and  poor  onsets  were  discarded  from  the  dataset,  along  with  those  events  too  poorly 
recorded  for  an  adequate  station  distribution  in  the  relocation.  That  the  average  absolute  rms 
residuals  drop  from  ~1.25  to  ~0.53  sec  Is  a measure  of  the  success  of  this  technique,  and  sug- 
gests that  about  half  the  uncertainty  in  the  absolute  locations  In  this  region  derives  from  the 
assumed  earth  model.  Ten  other  events  were  found  on  this  fault  trace.  Including  one  of  14  Sep- 
tember 1968,  magnitude  5.8,  which  had  a mechanism  almost  Identical  to  the  master,  and  are 
included  in  Fig.  IIl-9(c).  These  results  strongly  indicate  that  the  NE  dipping  plane  of  the  focal 
mechanism  is  the  fault  plane  and  that  the  relative  depths  obtained  are  meaningful.  It  Is  Inter- 
esting to  note  that  the  main  shock  is  at  the  top  of  the  sequence  which  extends  downward  some 
40  km.  Vast  salt  deposits  are  known  In  this  part  of  the  Zagros,  perhaps  preventing,  with  their 
ductility,  any  rupture  from  reaching  the  surface  (no  surface  break  was  found).  Another  ap- 
parent anomaly  is  that  the  dip  of  the  aftershock  zone  (60*  to  70*)  is  considerably  steeper  than 
that  of  the  fault-plane  solution  (36*)  with  a focal  velocity  of  6.8  km/sec.  The  dip  of  this  fault 
plane  in  the  focal  mechanism  increases  as  the  focal  velocity  decreases,  but  even  at  Vp  * 

6.0  km/sec  (the  minimum  likely  at  shallow  depths)  the  dip  increases  to  only  43*.  If  the  depths 
found  here  are  real,  it  would  suggest  that  either  the  fault  Is  steeper  at  depth  than  at  the  Initial 
break,  or  that  the  main  shock  oocurred  In  a limited  region  of  anomalously  low  velocity. 

The  effectiveness  of  this  technique  with  crustal  earthquakes  lends  it  to  several  applications. 
As  well  as  the  obvious  possibilities  of  resolving  fault-plane  ambiguities  In  focal  mechanisms, 
fault -plane  areas  may  be  estimated  and  are,  in  general,  better  defined  than  Is  possible  with 
either  geographic  or  conventional  joint  epioenter  determinations  which  use  residuals  in  the 
master-event  location  as  station  corrections  in  relocating  aftershocks.  If  aftershock  sones 
are  more  condensed  and  shallower  than  conventional  relocations  show,  the  corresponding  stress 
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Uiops  and  seismic  displacements  increase,  decreasing  the  discrepancy  between  seismic  and 

Q Q 

plate  motions  in  AbIo  * though  probably  not  removing  it  entirely.  This  work  has  been  contin- 
ued with  other  Iranian  aftershock  sequences,  and  may  he  extended  to  Central  Asia  in  the  future. 

T.J.  Fitch 
J.  Jackson* 


D.  TWO-DIM ENSIONAl.lZKD  VFI.OCITY  ANOMALIES 

1 ateral  heterogeneity  in  the  crust  and  upper  mantle  limits  the  accuracy  of  earthquake  lo- 
cations and  magnitude  determinations.  It  is  common  practice  to  add  station  corrections  to 
standard  travel  times  to  compensate  for  disparities  between  the  real  transmission  paths  and 
paths  inferred  from  a standard  earth  model,  thereby  achieving  more  accurate  locations.  Sim- 
ilarly, estimates  of  Mg  and  m^  can  be  corrected  for  station  amplitude  biases.  Both  corrections 
are  statistical  in  nature  and  thus  ignore  tile  underlying  physicB  of  seismic  transmission  in  the 
real  earth. 

Here  wo  present  a scheme  for  the  conversion  of  differential  arrival  times  to  velocity  anom- 
aly patterns  pertinent  to  near-source  regions.  With  teleseismic  data,  only  the  velocity  anom- 
alies seen  by  diving  rays  arc  resolved.  In  principle,  the  anomaly  patterns  can  be  used  to  com- 
pute near-source  corrections  to  magnitude  and  moment  as  well  as  to  travel  times. 

The  analysis  begins  by  locating  a group  of  events  relative  to  a master.7*10  Residuals  at 
each  station  are  averaged,  and  the  mean  residuals  are  projected  onto  a vertical  plane  normal 
to  the  trend  of  the  seismic  z-,ne.  This  is  carried  out  by  projecting  the  tangents  to  the  ray  paths 
from  the  master  into  the  appropriate  vertical  plane.  Angles  of  incidence  were  computed  from 
Herrin  and  J-B  earth  models  for  P and  PKP  phases,  respectively.  The  azimuths  of  the  ray 
paths  are  taken  from  absolute  locations. 

The  velocity  anomalies  are  linearly  relatod  to  the  mean  projected  residuals  through  velocity 
partial  derivatives.  The  assumption  is  made  that  classical  ray  parameters  computed  from  the 
assumed  earth  models  are,  at  least  on  the  average,  unchanged  by  lateral  heterogeneity  in  the 
source  region.  In  the  results  presented  here,  the  relative  locations  are  not  adjusted  for  the 
nonuniformity  of  the  source  velocity.  Consequently,  the  anomaly  patterns  may  contain  signifi- 
cant bias. 

Datasets  were  compiled  from  P and  PKP  arrivals  listed  in  the  ISC  Bulletins.  Records 
from  WWSSN  stations  in  the  region  of  the  activity  were  reread  if  the  reported  times  yielded 
travel-time  residuals  greater  than  4 to  5 sec. 

Anomaly  patterns  for  compressional  velocity  in  the  near-downward  vertical  direction  are 
presented  for  two  groups  of  earthquakes:  deep  focus  activity  within  130  km  of  a master  at  a 
depth  of  650  km  in  the  Tonga  seismic  zone  IFlgs.  Ill- 1 0(a)  and  (b)),  and  intermediate  depth  ac- 
tivity between  38*  and  40*N  beneath  Northern  Honshu  [Figs.  III-ll(a)  and  (b)].  Figure  III- 1 0(b) 
shows  a broad  velocity  anomaly  of  about  0.3  km/sec  in  the  down-dip  direction  of  the  seismic 
zone.  The  anomaly  pattern  for  the  Northern  Honshu  earthquakes  (Fig. III-ll(b)j  shows  a sharp 
minimum  and  maximum,  each  with  a magnitude  of  about  0.5  km/sec.  The  minimum  is  nearer 
the  down-dip  direction.  For  this  region  no  ray  paths  to  teleseismic  stations  leave  the  source 
region  in  the  down-dip  direction,  about  36*  toward  the  west. 

The  statistical  significance  of  the  anomaly  patterns  in  Figs.  111-10  and  in-11  is  tested  in 
Figs.  111-12  and  111-13.  By  definition,  a least-squares  solution  yields  residuals  with  a zero 

t Cambridge  University,  Cambridge,  England, 
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mean.  Con.si-quctitly,  the  mean  ol  any  subset  of  these  residuals  will  tend  to  zero  as  the  size  of 
the  subset  Increases,  If  the  residuals  are  random  variables.  This  trend  is  represented  by  the 
curves  in  I-'lga.  Ill- 1 2 and  111-13.  Data  points  above  these  curves  have  a significant  nonrandom 
component.  In  Fig.III-12,  mean  projected  residuals  from  the  Tonga  dataset  are  compared  with 
residuals  from  a synthetic  dataset  which  has  the  gross  statistical  properties  of  the  real  data. 

In  addition,  the  synthetic  data  were  computed  from  station  and  source  distributions  that  are 
Identical  to  the  real  ones.  Clearly,  the  synthetic  residuals  show  a higher  degree  of  random- 
ness; however,  several  synthetic  data  points  fail  the  test  for  randomness,  probably  because  of 
nonuniformity  in  the  spatial  distribution  of  sources  and  receivers. 

The  two-dimensionality  of  the  Honshu  seismic  zone  is  tested  in  Fig.  111-13.  The  data  points 
that  lie  above  the  2 a curve  correspond  to  the  maximum  and  minimum  in  the  anomaly  pattern  in 
Fig.  Ill- 1 1 . The  residuals  are  almost  completely  randomized  by  projecting  the  mean  residuals 
onto  a plane  45*  to  the  trend  of  the  zone.  However,  to  some  extent  this  effect  may  be  encouraged 
by  the  teleselsmlc  station  distribution  which  is  grossly  nonuniform  in  azimuthal  coverage. 

More  datasets  will  have  to  be  processed  to  establish  features  of  these  anomaly  patterns  that 
can  be  considered  characteristic  of  mantle  seismic  zones  in  different  depth  ranges  and  crustal 
activity  of  various  kinds.  Once  the  credibility  of  these  patterns  is  established,  then  amplitude 
effects  can  be  modeled.  j.  j 

J.  Jackson* 


E.  PRECISE  LOCATION  OF  OCEANIC  RIDGE  EARTHQUAKES 

The  world's  oceanic  ridge  system  is  characterized  by  Bhallow-depth  «5-  to  10-km)  earth- 
quakes of  intermediate  magnitude  (mb  4 to  6).  Detailed  knowledge  of  the  location  of  epicenters 
relative  to  eAch  other  and  to  the  ridge's  islands,  bathymetric  features,  and  magnetic/gravity 
fields  would  be  useful  for  geophysical  studies  as  well  as  nuclear  monitoring  activities.  Un- 
fortunately, conventional  methods  for  determining  earthquake  epicenter  locations  may  not  be 
adequate  to  attain  the  required  precision  and  accuracy,  in  that  relatively  simple  Jeffreys-Bullen11 
or  Herrin12  earth -velocity  models  are  used  to  calculate  eplcentral  distances.  These  models, 
which  generally  show  a thick  low  velocity  (6.0  km/sec)  continental  crust  overlying  a higher 
velocity  (7.8  km/Bec)  mantle,  are  clearly  not  applicable  to  oceanic  ridges.  Here,  high-velocity 

mantle  material  probably  extends  to  within  3 to  5 km  of  the  seafloor.  Also,  significant  velocity 

13  14 

inversion  may  exist  at  depth.  * 

The  computation  of  the  location  of  earthquakes  which  are  recorded  by  many  observatories 
distributed  in  an  azlmuthally  uniform  pattern  around  the  epicenter  is,  of  course,  insensitive  to 
the  errors  of  any  radially  symmetric  earth-velocity  model,  in  that  the  distance  error  along 
one  ray-path  azimuth  is  compensated  by  equal  error  in  the  opposite  direction  along  the  recip- 
rocal ray-path  azimuth.  By  simply  averaging  many  ray  paths,  a good  estimate  of  the  epicenter 
location  is  possible.  Unfortunately,  few  earthquakes  are  ever  recorded  by  such  an  Ideal  station 
distribution.  Thus,  epicenter -location  errors  can  be  expected  to  be  biased  in  the  general  direc- 
tion toward  the  greatest  number  of  recording  stations. 

This  biasing  effect  is  well  exhibited  by  the  earthquake  epicenters  which  have  been  located 
along  the  central  mid-Atlantic  ridge  by  the  ISC  for  the  years  1964  to  1973  and  the  USGS  since 
1973.  Here,  the  relatively  small  ridge-crest  events  are  primarily  recorded  by  North  American 
observatories  (n  ~ SO),  since  there  are  few  North  African  observatories  (n  **  5). 

t Cambridge  University,  Cambridge,  England. 


48 


I igtii  «■  111-14  shows  Hu-  ISC'  l umputod  epicenter  locations  relative  t.i  the  topographically 
defined  motlian-i  lft  valley  and  the  transform  fault  tectonic  axis  between  15*  and  29*N  latitudes. 

The  topography  information  was  taken  from  recently  compiled  detailed  U.S.  Navy  bathymetry 
contour  charts  (Naval  Research  l.ab< ..  «iorv).  The  positional  accuracy  of  the  axial  tectonic  loci 
is  estimated  to  be  on  the  order  of  *4  km.  For  this  study,  the  tectonic  axis,  as  defined  by  the 
bathymetric  data,  has  been  digitized  for  incorporation  into  a variety  of  plotting  and  analysis 
programs.  The  digitization  process  was  applied  over  small  individual  segments  of  the  ridge, 
whose  si/e  dejiended  on  the  detail  of  the  bathymetric  charts  available.  When  completed,  the 
mid-Atlantic  ridge  axis  will  have  been  continuously  digitized  between  5*  and  65*N  latitudes.  The 
total  digital  length  of  the  ridge  is  expected  to  produce  a data  base  of  over  30,000  points,  or  about 
2 points  each  kilometer  of  median  valley  and  transform  fault  length.  While  this  degree  of  resolu- 
tion may  not  be  necessary  in  the  initial  stages  of  the  study  reported  here,  It  will  be  required 
later  for  developing  detailed,  local  ridge  axis,  crustal/upper-mantle  models. 

Careful  inspection  of  Fig.IlI-14  reveals  that  the  epicenters  are  rather  widely  scattered,  but 
most  arc  clearly  clustered  west  of  the  northerly  trending  median  valley  segments  or  north  of 
the  easterly  trending  transform  faults.  One  might  have  expected  the  epicenters  to  have  been 
normally  distributed  about  the  tectonically  active  median  valley  and  transform  fault  axes.  In 
fact,  detailed  microoarthquake  studies  on  the  mid-Atlantic  rtdgc  islng  local  ocean-bottom  seis- 
mometers and  sonobuoys  show  all  seismicity  to  be  confined  to  within  5 km  of  the  median  valley 
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and  transform  fault  axes.  This  systematic  location  offset  dlscre;>ancy  of  about  15  to  20  km, 

which  is  best  seen  in  the  vicinity  of  the  Kane  Transform  Fault  at  24 *N  latitude  (Fig.  III-15), 
probably  reflects  the  simple  fact  ‘hat  most  of  the  recording  stations  are  Indeed  located  to  the 
north  and  west  of  the  ridge  in  North  America,  Hius,  use  of  the  relatively  "slow"  J-B  earth- 
velocity  model  by  the  ISC  rather  than  a more-appropriate  "fast"  oceanic  ridge  axla  model  would 
be  expected  to  "pull"  the  epicenters  toward  the  biasing  North  American  stations.  An  Increase 
in  the  average  cr  jatal  velocity  beneath  the  ridge  axla  by  about  20  percent  la  Implied. 

The  large  scatter  in  the  epicenter  locations,  especially  along  the  transform  faults,  as  com- 
pared with  the  narrow  tectonic  activity  zone  inferred  by  the  local  microearthquake  studies.  Is 
alBo  noteworthy.  The  scatter  may  be  real,  which  would  suggest  a fundamental  difference  between 
small  microearthquake  and  large,  teleselsmically  recorded  earthquake  mechanisms.  This  Is 
undoubtedly  true  for  some  of  the  large  teleseismlc  events  located  well  off  the  ridge  axis  O 50  to 
100  km)  which  show  thrust  mechanisms.  Alternatively,  the  scatter  could  simply  result  from 
the  same  station-biasing  effect  which  has  offset  the  entire  epicenter  pattern  from  the  ridge  axis, 
t may  be  that  relocation  of  the  ISC  and  USCS  epicenters  with  an  improved  "oceanic  ridge"  ve- 
locity model,  which  translates  the  epicenter  pattern  onto  the  median  valley  and  transform  fault 
axes,  will  also  reduce  the  noar-axis  scatter. 

J.  D.  Phillips 
R. M.  Sheppard 
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Fig.  111-1.  Result  of  convolving  least-squares 
Inverse  filter  with  SRO  SP  impulse  response. 
Filter  is  2.0  sec  long  and  was  designed  with 
5-percent  additive  white  noise.  Output  occurs 
with  delay  of  1.4  sec. 
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Fig.  111-2.  Results  of  deconvolving  and  anti-echo  filtering  MAIO  recording 
of  20  March  1076  E,  Kazakh  event. 
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Fig.  111*3.  Maximum  entropy  cepatrum 
of  deconvolved  seismogram. 
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Fig.  III-8#  Plot  of  reported  depth  (ISC)  against  number  of  stations 
recording  event.  Note  that  few  deeper  events  are  well  recorded. 
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Ftg.UI-9.  (a,  b)  Relocation  of  1972  QIR  sequence.  Error  bar* 
are  standard  errors  checked  by  xz  tests.  Master  event,  labeled 
M,  la  at  28.39  “N,  52.78*E,  11  kin.  Those  events  with  poor  sta- 
tion distributions  are  labeled  □ . (c)  Equal-area  plot  of  position 
vectors  to  aftershocks  on  focal  net  of  master  event.  Compare 
trend  with  mechanism  in  (a).  This  is  plot  of  all  located  activity 
on  fault  trend. 
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Kig.  Ill- 14.  Earthquake  distribution  along  a portion 
of  mid-Atlantic  ridge. 


IV.  DATA  ANALYSIS  AND  COMPUTER  SYSTEMS 


A.  TRANSFER  FUNCTIONS  FOR  TIIE  SEISMIC  RESEARCH  OBSERVATORY 

SEISMOGRAPH  SYSTEM 

Tlie  Seismic  Research  Observatories  (SROs)  can  produce  long-period  (LP),  short- 
period  (SP).  or  broadband  3 -component  data  which  can  be  either  analog  or  digital.  It  is  likely 
that  the  data,  particularly  the  digital  data,  will  be  widely  used  for  seismic  research;  therefore, 
it  is  important  that  response  ■ haracterlstics  be  well  known  and  understood.  This  is  particularly 

important  because  these  seismographs  are  quite  different  from  classical  seismic  instrument 

. 1.2 

systems. 

It  is  our  objective  to  present  here  all  important  dynamic  systems'  response  information 

which  may  be  needed  to  understand  or  interpret  SHO  data.  Our  approach  is  basically  that  used 

to  design  the  system.  We  derive  our  results,  as  the  designing  engineer  did,  from  either  circuit 

3 4 

diagrams  or  mathematical  specifications.  While  more  involved  than  approximate  methods,  ' 
this  procedure  allows  the  overall  system  response  to  be  divided  into  individual  components. 

No  distinction  is  made  between  vertical  and  horizontal  instruments  because  they  have  the  game 
response  functions.  * ^ 

The  models  we  have  constructed  are  based  on  information  in  SRO  manuals.  In  addition,  a 
number  of  individuals  and  companies  have  been  very  helpful  to  us  in  obtaining  the  SRO  informa- 
tion we  needed.  We  particularly  want  to  thank:  the  White  Instrument  Company  of  Houston,  Texas 
for  information  about  the  anti-alias  filters;  the  Teledyne-Geotech  Company  of  Garland,  Texas 
for  details  of  the  seismometer  and  shaping  filters;  Mr.  Howard  Durham  of  Sandia  Laboratories, 
Albuquerque,  New  Mexico  for  additional  inputs  covering  the  seismometer;  and  Mr.  Gary  Holcomb 
of  the  USGS  Albuquerque  Seismological  Center  who  has  been  very  cooperative  and  helpful  in  ob- 
taining information  as  we  required  it. 

It  le  important  to  note  that  our  models  are  general,  but  the  SRO  SP  systems  now  operate 
without  an  anti-alias  filter.  It  appears  that  the  anti-aliaB  filter  was  removed  at  all  operating 
sites  during  the  first  half  of  1976,  but  at  different  times.  The  correct  SRO  SP  response  may  or 
may  not  include  the  anti -alias  filter,  depending  on  what  date  the  recording  was  made. 


i.  The  Seismometer 


The  basic  SRO  seismometer  is  a Teledyne -Geotech  Model  36000  which  includes  a very 
lightly  damped  accelerometer  with  a natural  frequency  of  approximately  0.2  Hz.  A capacitance 
displacement  measurement  generates  a feedback  signal  which  is,  in  turn,  used  to  maintain  null 
position.  This  instrument  Is  therefore  similar  In  design  to  the  lunar  seismometer  analyzed  by 
Sutton  and  Latham.5 

The  seismometer  has  two  basic  outputs.  One,  v^,  is  a broadband  output  which  is  essen- 
tially flat  to  acceleration  inputs  from  DC  to  1.0  Hz.  Above  1.0  Hz,  the  response  drops  at  12  dB/ 
octave.  The  other,  vf,  Is  a filtered  response  which  Is  similar  except  that  additional  filtering 
causes  the  response  to  drop  at  6 dB/octave  at  frequencies  below  0.02  Hz.  Mere  exactly,  the 
acceleration  to  broadband  output  transfer  function  is 


vbb(9) 

-tfe-  * Hbba<9> 


-197  (s  + 0.125)  (s  + 50) 
(TTOW1S  +0277T+tent  + 3 57IT 


(IV-l) 


and  the  filtered  transfer  func  tion  in 


vf(8)  ..  , , 39,400s  (s  + 50) 

Tun  n ' fa  * = (7T0.tWT^+  IM\U  ♦ UM  (a  ♦ tt.f) 


(IV -2) 


where  output  Is  measured  in  volts,  and  acceleration  is  in  cm/sec^.  In  both  cases,  the  transfer 
function  for  displacement  input  is  obtained  by  multiplying  the  above  by  s*. 

Figure  IV-i(a-b)  shows  the  broadband  system  response  to  acceleration  inputs.  As  noted 
above,  the  amplitude  response  is  flat  from  DC  to  approximately  1.0  Hz.  The  phase  delay  is 
shown  assuming  that  the  1 in  corresponds  to  a phase  shift  of  +*  radians.  Figure  IV-2(a-b) 
shows  the  corresponding  filtered  response  to  displacement  inputs. 


2.  Response -Shaping  Filters 

The  SRO  system  produces  both  LP  and  SP  outputs  by  furthor  filtering  of  the  filtered  output 
Vf  of  the  basic  seismometer. 

a.  LP  Response  Shaping 

The  LP  shaping  filters  have  been  selected  to  yield  a system  response  to  earth  displacement 
having  a peak  response  at  a 28-aec  period,  a 1.25-dB  octave  separation  of  -3-dB  points,  a 
24-dB/octave  roll-off  at  high  and  low  frequencies,  and  a 6-sec  notch  to  remove  microseisms. 
This  has  been  done  by  cascading  the  first  four  filters  with  the  transfer  functions  shown  in 
Table  IV -1.  The  fifth  filter  is  also  in  the  system  to  block  DC  offsets;  it,  in  fact,  increases 
the  roll-off  to  30  dB/octave  at  very  long  periods.  The  overall  shaping  filter  transfer  function 
is  Just  the  product  of  the  transfer  functions  in  Table  1V-1. 


b.  SP  Response  Shaping 

The  SF  shaping  alters  have  been  selected  to  yield  a system  response  to  earth  displacement 
having  a peak  response  at  a 0.35-sec  period  (2.86  Hz),  a 1.8 -dB/octave  separation  of  -3-dB 
points,  and  an  18-dB  roll-off  toward  lower  frequencies  in  the  microseism  band.  The  filters 
applied  to  the  basic  seismometer  filtered  output  to  achieve  this  are  given  in  Table  IV-2.  Note 
that  for  frequencies  much  below  10  Hz,  the  10-Hz  bandpass  filter  might  as  well  be  a simple 
differentiator  giving  6 -dB/octave  loss  toward  lower  frequencies. 


3.  Anti -alias  Filters 

Anti-alias  tiltora  were  originally  applied  to  both  LP  and  SP  analog  outputs  prior  to  sampling 
and  digital  recording.  These  filters  are  introduced  to  be  certain  that  there  is  no  significant 
energy  remaining  above  the  Nyqulst  frequency.  The  anti -alias  filters  further  modify  the  instru- 
ment response,  particularly  close  to  and  above  the  Nyquist  frequency.  As  noted  above,  the  SP 
system  now  operates  with  no  anti -alias  filter  and  depends  upon  nature  to  produce  a signal  with 
no  energy  above  10  Hz  (sampling  rate  is  20  Ha). 


a.  LP  Anti-alias 

The  LP  anti -alias  Alter  is  a 4 -pole  Buvterworth  low -pass  Alter.  The  response  varies 
smoothly  from  0 dB  at  DC  to  -3  dB  at  0.1  Hs,  and  then  drops  to  -60  dB  at  the  0.5-H*  Nyquist 
frequency. 
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The  form  of  the  filter  transfer  function  is 


IT  (a  + aw) 
kai  K 


(IV-3) 


where  u>0  Is  the  natural  (3-dB)  frequency  of  the  filter,  and  the  - • ^ are  the  poles  of  the  filter. 
The  values  for  the  LP  anti -alias  filter  constants  are: 


w0  ■ 0.1(2t ) > 0.628 


st  > 0.240  >0.5801 
»2  * 0.580  - 0.2401 
Sj  ■ 0.580  4 0.2401 
S4  ’ 0.240  4 0.5801 


b.  SP  Anti -alias 

The  original  SP  and-aliaa  filter  was  a 7-pole  Tchebyehev  low-pass  filter.  It  had  0.25-dB 
ripple,  was  down  3 dB  at  5.0  He,  and  down  more  than  60  dB  at  10.0  Hz.  The  small  ripple  and 
larger  number  of  poles  compared  with  the  LP  anti-alias  filter  were  required  because  of  the 
specification  of  almost  60 -dB  drop  in  the  single  octave  from  5.0  to  10.0  Hz. 

The  form  of  the  filter  was 

An  zkak i 

Ht  a Jyfat. (IV *4) 

(S  4|.| 
kal  * 

where  the  -sk  are  the  poles  of  the  filter.  The  numerical  values  for  these  were: 

St  > 2.055  -30.631 
S2  ■ 5.759  -24.561 
a,  > 8.321  - 13.631 
s4  ■ 9.236  - 0.01 
s4  * 8.321  4 13.631 
l6  ■ 5.759  4 24.561 
s?  « 2.055  4 30.631 

The  resulting  numerator  constant  In  H(1  selected  to  give  unity  gain  at  DC,  was  1.413  x 10^. 

4.  Overall  Response  Characteristics 

The  transfer  functions  in  the  previous  section  are  complete  linear  descriptions  of  the  LP 
and  SP  responses  of  the  SRO  seismic  system.  We  can  now  consider  the  overall  response  char- 
acteristics. Figure  IV-3(a-b)  shows  the  SRO  LP  response  to  displacement,  and  Flg.lV-4(a-b) 


Allow*  the  SI*  response  to  displacement  with  and  without  the  anti-alias  filter*.  In  both  case*,  it 
in  clear  that  the  anti-alias  filters  nitniflcantly  modify  the  high-frequency  amplitude  reeponeee 
and  that  there  are  sul>ntuntial  phase  > ontributions  well  Into  the  usual  frequency  banda  of  Interest. 
Figures  IV*<i  and  IV-H(a-e)  show  the  corresponding  system  Impulse  responses. 

An  important  characteristic  of  seismograph  systems  is  the  group  delay  they  introduce. 

This  is  the  negative  of  the  derivative  of  the  transfer  function  phase: 


y-> » 


(IV—  5) 


where  o(w)  is  the  phase  response.  We  have  calculated  t^  for  the  SRO  LP  transfer  function  by 
numerically  differentiating  4(u>).  The  resulting  group  delay  curve,  computed  with  and  without 
the  anti-alias  filter  for  purposes  of  comparison,  is  shown  in  Fig.  IV-5.  Since  the  LP  group 
delay  is  of  interest  for  a number  of  applications,  we  have  derived  a simple  parametric  repre- 
sentation.  This  approximation,  which  Includes  the  effect  of  the  anti-alias  filter,  is 

t 

(I.M1  K iO‘^1  ♦ (t  .96(1  x tO**u|  - 0.166  X ♦ (9.9SS  x - (*.511  x I0'*u>4) 
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It  gives  a root-mean-square  error  of  less  than  0.5  sec  between  a 10-  and  150-sec  period.  The 
worst  fits  are  at  the  longer  periods,  and  are  as  large  as  1.25  sec. 

Part  of  the  group  delay  of  the  LP  instrument  is  due  to  the  anti-alias  filter.  For  example, 
there  is  an  apparent  pure  delay  of  a few  seconds  on  the  impulse  response  shown  in  Fig.  IV -6. 
That  is  primarily  due  to  the  anti-alias  niter.  In  general,  the  LP  anti-alias  filter  does  little 
more  than  introduce  that  delay  and  slightly  smooth  the  impulse  response.  We  will  see  below 
that  this  is  not  the  case  for  the  SP  system. 

It  is  not  customary  to  calculate  or  discuss  group-delay  curves  for  SP  instruments.  How- 
ever, the  SRO  SP  is  special  and  worth  considering.  Figure  IV-7  shows  the  SP  group  delay  for 
the  SRO  and  WWSSN  SP  instruments.  The  transfer  function  for  the  WWSSN  is  given  by 
Chakrabarty  fiii2  Except  near  1 Hz,  the  SRO  with  an  anti -alias  filter  clearly  introduces 
more  delay.  Moreover,  the  SRO  actually  has  a group  delay  minimum  which  corresponds  to  a 
buildup  of  3-Hz  energy  with  a delay  of  some  0.35  sec.  The  WWSSN  has  small  high-frequency 
delays  giving  a sharp  onset,  while  the  SRO  with  anti -alias  filters  has  large  high-frequency 
delays  giving  an  emergent  onset.  The  impulse  responses  shown  in  Fig.  lV-8(a)  (SRO  with  anti- 
alias  filter),  IV-8(b)  (WWSSN),  and  IV-8(c)  (SRO  without  anti-alias  filter)  can  be  interpreted  in 
terms  of  the  group-delay  curves.  In  particular,  the  maximum  at  0.35  sec  and  the  5-Hz  ringing 
which  comes  in  later  when  the  system  includes  an  anti-alias  filter  are  clearly  visible.  Users 
of  very  early  SP  SRO  digital  data  should  be  aware  of  these  artifacts  in  the  Instrument  response 
in  order  to  properly  interpret  their  data. 

It  should  be  noted  that  the  current  configuration  (no  anti-alias  filter)  for  the  SRO  SP  system 
can  cause  some  difficulties.  Background  noise  and  most  teleseismic  signals  have  relatively 
little  energy  at  frequencies  above  a few  Herts.  There  is  usually  little  energy  above  10  Hs  to  be 
aliased,  so  most  teleseisms  should  be  faithfully  recorded.  However,  the  instrument  response 
is  not  down  very  much  at  10  Hz,  so  local  noise  or  events  with  much  more  high-frequency  content 
might  experience  aliasing.  It  might  be  prudent  to  introduce  some  additional  filtering  into  the 
SP  system,  but  not  a filter  as  extreme  as  the  original  anti -alias  filter. 


R,  T.  Lacoss 
D.W.  Me  Cowan 


H.  DETECTION  RESULTS  l'OR  AN  SRO  NETWORK 

The  development  and  installation  of  a Seismic  Research  Observatory  Network  has  proceeded 
to  a point  that  now  allows  for  some  analysis  of  the  event  detection  performance  of  the  observe > 
tories.  In  the  normal  dally  operations  of  the  SROs,  an  event  detector  is  attached  to  the  output 
of  the  seismometer  SP  vertical  channel.  Events  which  are  sufficiently  large  to  activate  the 
detector  will  cause  the  SP  output  of  the  seismometer  to  be  recorded  on  magnetic  tape.  A com- 
plete list  of  the  event  detections  and  recorded  time  periods  is  then  maintained  in  a detection  log. 
About  every  two  weeks,  the  data  tapes  are  shipped  to  the  USGS.  Once  the  SRO  data  tapes  arrive 
at  Albuquerque  Seismological  Laboratory,  they  are  processed  and  checked  for  data  integrity. 

The  detection  logs  and  the  actual  SRO  digital  data  are  then  transmitted  over  the  ARPANET  to 
the  mass  store  Datacomputer.  Within  the  Datacomputer,  the  detection  logs  are  stored  in 
monthly  file  summaries  while  the  digital  data  are  stored  in  individual  day  files.  Also  within 
the  Datacomputer  are  monthly  bulletin  files  which  have  been  compiled  from  stations  that  rapidly 
report  their  arrival  times  to  the  USGS.  The  bulletin  files,  called  PESF  or  Preliminary  Event 
Summary  Files,  are  meant  to  provide  a very  preliminary  location  only  to  indicate  the  time 
periods  of  SRO  data  that  may  be  of  further  interest.  It  is  the  PESF  and  the  SRO  detection  log 
files  that  form  the  basis  of  this  study. 

A comparison  of  the  event  detections  of  an  SRO  network  of  stations  consisting  of  ANMO, 
CTAO,  GUMO.  MAJO,  NWAO,  and  ZOBO  and  the  event  bulletins  in  the  PESFs  of  the  Datacom- 
puter reveal  some  interesting  figures.  There  were  only  six  SROs  in  operation  during  the  58-day 
study  period,  1 April  to  28  May  1977.  The  detection  logs  for  the  SROs  were  obtained  from  the 
detection  log  files  In  the  Datacomputer.  Event  locations  for  the  corresponding  period  were  also 
extracted  from  the  Datacomputer  PESFs.  A summary  of  the  monthly  event  totals  and  SRO  de- 
tections is  given  in  Table  IV-3  where  the  basic  data  as  extracted  from  the  Datacomputer  files 
are  collected  and  tabulated.  During  this  time  period,  most  of  the  SROs  operated  continuously, 
with  two  exceptions.  The  site  GUMO  was  inoperative  for  the  last  half  of  April  and  the  first  5 
days  of  May.  CTAO  was  inoperative  during  the  entire  month  of  May.  For  the  remaining  SROs 
there  is  a wide  variation  in  the  number  of  event  detections  reported  during  the  study  period. 

For  example,  MAIO  reported  427  detections,  while  ZOBO  reported  918  for  the  same  time  pe- 
riod. NWAO  reported  374,  while  ZOBO  reported  1138  for  the  same  28-day  period  in  May. 
Clearly,  the  event  detection  thresholds  are  not  in  agreement  between  various  SROs.  Even  ac- 
counting for  differences  in  background  noise,  the  sites  are  not  recording  detections  at  consistent 
detection  levels. 


TABLE  IV-3 

NUMBER  OF  SRO  DETECTIONS 

Month 

Bulletin 

Event* 

ANMO 

CTAO 

OUMO 

MAIO 

NWAO 

ZOBO 

April 

326 

332 

437 

BO 

427 

176 

918 

Moy  (1-28) 

287 

632 

0 

993 

394 

374 

1138 

TABLE  IV-4 

EVENTS  DETECTED  AS  A PERCENT  OF  TOTAL  BULLETIN  EVENTS 

Month 

ANMO 

CTAO 

GUMO 

MAIO 

NWAO 

ZOBO 

April 

52% 

30% 

46% 

45% 

16% 

50% 

Moy  (1-28) 

56% 

0% 

44% 

45% 

21% 

63% 

During  the  58-day  study  period,  the  performance  of  the  SROb  varied  widely.  Using  the 
bulletin  event  parameters  of  location  and  depth,  the  arrival  tim>‘j  were  computed  for  each  of 
the  SRO  sites.  If  an  SRO  detection  fell  within  the  computed  event  arrival  time  *4  min.,  the 
event  was  considered  bb  being  detected  at  that  SRO.  In  Table  IV -4  the  events  detected  are  ex- 
pressed as  a percent  of  the  total  events  in  the  bulletin  and  are  listed  for  each  SRO  for  both 
months.  From  this  table,  it  is  clear  that  there  are  two  good  and  two  moderately  good  SROs, 
and  one  poor  SRO.  Because  of  the  extensive  down  time,  data  for  CTAO  are  insufficient  to  draw 
any  conclusions.  Both  ANMO  and  ZOBO  are  operating  with  a good  level  of  detections,  and  are 
recording  over  50  percent  of  the  bulletin  events.  QUMO,  when  operating,  performs  as  well  as 
MAIO,  both  detecting  somewhat  less  than  50  percent  of  the  bulletin  events.  The  one  site  that 
is  in  need  of  improvement  is  NWAO,  although  increasing  the  number  of  detections  at  NWAO  by 
lowering  the  detection  threshold  may  not  result  in  any  great  increase  in  real  event  detections. 
For  example,  comparing  NWAO  and  MAIO  during  May  indicates  that  both  sites  had  about  the 
same  number  of  total  detections,  yet  for  NWAO  the  real  events  (bulletin  events)  detected  are 
about  half  those  of  MAIO. 

Based  on  the  total  number  of  detections  listed  for  these  two  months,  it  appears  that  moBt 
of  the  SROs  are  operating  with  a detection  threshold  that  is  set  too  high.  Lowering  the  threshold 
at  NWAO,  CTAO,  GUMO,  and  MAIO  would  increase  the  total  events  detected  by  the  network. 

A clear  example  of  this  is  seen  in  the  performance  of  ZOBO. 

R.  M.  Sheppard 


C.  AUTOMATIC  PICKING  OF  FIRST  ARRIVALS 

A simple  algorithm  is  being  studied  which  Identifies  earthquakes  and  times  first  arrivals 
using  data  from  a local  network  of  single -component  seismometers.  The  algorithm  is  simple 
enough  that  it  could  easily  be  Implemented  on  a microprocessor  monitoring  a single  station  in 
real  time. 

A unique  aspect  of  the  algorithm  is  that  it  not  only  picks  arrivals,  but  also  objectively 
determines  arrival-time  accuracy  after  arrivals  have  been  picked.  Earthquakes  are  located 
by  an  Iterative  method.  First,  the  most -accurate  arrivals  are  used  to  locate  the  event.  Then, 
arrivals  which  do  not  agree  with  this  location  are  replcked  using  more -sensitive  parameter 
settings.  This  process  is  repeated  until  a stable  location  is  obtained.  Using  such  an  algorithm 
in  routine  processing  of  seismic  data  could  require  as  little  as  2 percent  of  the  human  effort 
now  required  by  manual  processing  and  could  provide,  for  the  first  time,  an  accurate,  objective 
data  base  for  earthquake  research. 

The  simplest  algorithm  that  seemed  to  work  reliably  in  most  cases  was  developed.  The 
idea  behind  this  was  that  the  analysis  of  bugs  in  this  algorithm  will  show  how  the  algorithm 
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could  be  improved.  In  fact,  the  bugs  thut  were  found  are  general  enough  that  they  must  be 
contended  with  by  any  picking  algorithm.^ 

There  is  always  an  ambiguity  associated  with  measuring  the  first-arrival  time.  This  is 
due  to  the  fact  that  the  signal  is  contaminated  with  noise,  the  shape  of  the  arrival  is  unknown, 
and  the  arrival  waveform  varies  from  station  to  station.  Thus,  the  methods  of  "beam  forming" 
and  "matched  filtering"  are  unfortunately  not  applicable. 

However,  the  first  arrival  does  occur  in  a context,  the  features  of  which  can  be  determined 
more  reliably.  These  features  can  be  used  to  suggest  where  the  first  arrival  should  be,  and 
also  Increase  the  reliability  of  the  pick.  The  algorithm  uses  a simple  first-arrival  detector 
followed  by  additional  semantic  checks  to  be  sura  that  the  pick  is  associated  with  a seismic 

event.  The  details  of  the  algorithm  are  presented  in  Ref.  6.  The  algorithms  of  Stewart,7 
8 9 

Stevenson,  and  Allen  operate  similarly. 

The  data  used  to  test  the  algorithm  are  from  a magnitude  3.5  earthquake  recorded  by  the 
Central  California  seismic  network;  47  seismic  tracers  were  used.  Examples  of  the  arrival 
times  determined  by  the  algorithm  are  shown  in  Fig.  IV  -9.  The  vertical  dotted  line  1b  the  com- 
puted arrival  time,  and  the  number  to  the  right  of  it  is  an  estimate  of  the  standard  error  of  the 
arrival  time  in  seconds.  For  each  trace,  1 sec  of  noise  before  and  1 sec  of  signal  after  the 
arrival  time  are  shown.  The  tick  marks  are  0.1  sec  apart.  The  portions  of  the  seismograms 
corresponding  to  large  amplitudes  have  been  removed  to  clarify  the  figure.  The  traces  were 
chosen  to  show  how  the  algorithm  performed  under  various  noise  and  signal  conditions.  Note 
the  large  variations  in  the  first-arrival  waveform. 

Of  the  47  seismic  signals,  3 were  correctly  classified  as  telephone -line  noise.  The  arrival 
times  for  sharp  arrivals  usually  agreed  to  within  a hundredth  of  a second  of  hand  picks.  For 
weaker,  more  ambiguous  arrivals,  computer  picks  and  hand  picks  are  always  within  the  esti- 
mated error  calculated  by  the  algorithm. 

The  programming  for  a complete  automated  system  has  not  been  finished,  but  the  analysis 
might  proceed  as  follows.  First,  the  picking  algorithm  would  be  operated  on  all  the  seismic 
channels.  If  the  spatial  and  temporal  distribution  of  arrivals  over  the  network  was  such  that 
an  event  was  likely  to  have  occurred,10  the  hypothesis  that  an  event  had  occurred  would  be 
pursued.  If  there  were  enough  accurate  arrivals,  an  attempt  would  be  made  to  locate  the  event. 
For  the  earthquake  studied,  17  arrivals  were  given  accuracies  of  at  least  0.01  sec,  and  these 
were  used  to  locate  the  event.  The  location  found  agreed  well  with  that  determined  from  hand- 
picked arrivals.  The  arrivals  would  then  be  re-examined  with  a more-sensitive  algorithm  to 
adjust  the  pick  if  necessary.  Arrival-time  residuals  could  provide  some  information  to  guide 
the  processing  at  this  point.  For  example,  a positive  residual  indicates  that  the  arrival  is  late 
compared  with  what  the  velocity  model  predicts.  Thus,  if  the  station  is  relatively  flar  from  the 
earthquake  (greater  than  70  km,  say),  this  suggests  that  the  algorithm  should  look  for  a weak 
longer-period  arrival  that  might  be  a refracted  arrival  originally  missed.  A large  negative 
residual  means  that  the  arrival  is  earlier  than  predicted,  and  might  suggest  that  a glitch  may 
precede  the  actual  arrival. 

However,  one  cannot  rely  on  residuals  too  heavily  since  they  can  reflect  more  the  inade- 
quacies of  the  model  than  inaccuracy  of  the  picks.  For  example,  although  10  arrivals  have 
residuals  greater  than  0.5  sec,  Station  RU8's  are  the  only  true  gross  errors.  RUS  is  the 
farthest  station  reporting  an  arrival,  and  has  a 2. 13 -sec  residual  suggesting  that  a refracted 
arrival  was  missed.  RUS  was  repicked  using  a lower  trigger  threshold,  and  a residual  of 
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0.13  sei  waa  obtained.  The  initial  and  final  picks  are  shown  in  Fig.  IV-10.  It  should  be  noted 
that  the  algorithm  dot’9  not  operate  reliably  using  such  a reduced  threshold. 

Thus,  the  algorithm  has  the  potential  of  using  partial  information  (in  this  case,  a tentative 
location)  to  operate  on  the  data  in  a way  that  could  not  be  used  otherwise.  Using  multiple  levels 
of  information  makes  processing  much  more  complicated.  However,  it  provides  constraints 
wtiich  are  not  available  at  the  lower  levels  of  the  system.  Using  simple  tests  ("triggers")  to 
suggest  possibilities  allows  partial  knowledge  to  be  used  to  control  the  analysis  in  a manner 
which  is  not  irrevocable.  That  is,  an  hypothesis  can  be  suggested,  and  either  accepted  or  re- 
jected, and  processing  modified  accordingly.  Thus,  the  algorithm  is  truly  data  adaptive. 

K.  R.  Anderson 


D.  PDP-11  SYSTEM  SOFTWARE 

During  the  last  six  months,  we  have  continued  to  develop  our  PDP-11  UNIX  software  system 
to  meet  the  needs  of  the  Group.  Three  important  developments  have  been  the  installation  of  a 
new  magnetic -tape  driver,  the  development  of  a new  software  maintenance  system,  and  initial 
work  on  Multi-Computer  UNIX. 

The  new  magnetic-tape  driver  removes  three  major  defects  from  the  standard  UNIX  driver: 
Inability  to  use  two  tapes  simultaneously,  inability  to  process  multiple  files  on  one  tape,  and 
inability  to  recover  data  from  tapes  with  noise  and  errors. 

PDP-11  software  maintenance  has  become  a difficulty,  because  a large  volume  of  software 
hurriedly  imported  to  get  the  system  going  did  not  have  adequate  source  code  or  documentation, 
and  because  insufficient  on-line  storage  resulted  in  poor  organization  of  our  rather  extensive 
program  source  code  and  computer-readable  manual  libraries.  We  are  now  halfway  through  a 
project  to  put  all  source  code  and  documentation  on-line,  and  to  be  sure  that  the  source  code 
and  documentation  correspond  to  the  running  system  programs.  Several  Bimple  maintenance 
programs  are  being  developed  to  make  it  possible  to  keep  the  system  in  a crisp,  well-maintained 
state  after  this  initial  shakedown  is  completed.  The  shakedown  and  maintenance  programs  are 
to  be  finished  by  the  end  of  1977.  This  project,  which  requires  a considerable  effort  on  our 
part,  is  also  important  to  others  - such  as  the  VELA  Selsmologlcal  Center,  to  whom  we  have 
promised  our  software. 

The  load  on  our  PDP-11/ 50  running  UNIX  has  become  quite  heavy,  and  it  la  clear  that  as 
more-and-more  selsmologlcal  algorithms  are  programmed  for  the  computer,  and  as  work  is 
further  transferred  from  our  now  antiquated  PDP-7's,  one  PDP-11  will  no  longer  be  adequate. 

We  are  in  the  process  of  revising  UNIX  to  run  on  two  computers  simultaneously,  so  that  we  can 
use  our  PDP-11/40  to  lighten  the  UNIX  load.  The  PDP-il/40  is  the  same  size  as  the  PDP-U/50, 
except  that  its  CPU  is  half  the  speed;  but  should  CPU  speed  become  a problem,  the  PDP-11/40 
can  be  easily  and  cheaply  upgraded  to  a PDP-U/60  which  is  slightly  faster  than  the  PDP-il/50. 
The  revised  system,  called  Multi-Computer  UNIX,  must  share  disks  and  teletypes  between  the 
computers.  Multi-Computer  UNIX  has  gone  through  Initial  design  and  about  half  of  initial  coding 
since  the  last  semiannual  report. 

R.  L.  Walton 

E.  UNIX  SIGNAL  DISPLAY  PACKAGE 

Work  is  being  done  to  implement  a general-purpose  interactive  seismic  Signal  Display 
Package  (SDP)  using  storage  scope  displays  to  run  on  our  PDP-11  UNIX  operating  system.  To 
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facilitate  this  work,  the  graphics  system  described  in  the  previous  SATS11  has  been  expanded 
to  include  a command  to  plot  a waveform  directly  from  a UNIX  file.  Originally,  the  data  would 
have  to  have  been  read  by  the  user  program  and  then  passed  to  the  graphics  routines.  This 
change  not  only  reduces  the  time  necessary  to  plot  the  waveforms,  it  also  saves  the  user  from 
having  to  allocate  space  for  the  large  data  arrays.11 

The  current  design  of  the  SDP  consists  of  four  parts:  (1)  the  seismic  waveform  database 
which  must  be  in  our  standard  UNIX  waveform  format;  (2)  the  general  display  or  interface  pro- 
gram; (3)  the  display  command  file  which  is  an  alphanumeric  file  of  applications -oriented  prim- 
itive graphics  commands  describing  the  display;  and  (4)  a graphics  command  Interpreter  which 
is  a subroutine  to  translate  the  primitive  commands  into  actual  graphic  commands. 

The  standard  format  for  UNIX  waveform  data  consists  of  an  alphanumeric  index  file  and  a 
group  of  binary  data  files.  The  Index  file  contains  one  line  of  formatted  Information  describing 
each  waveform,  and  the  data  files  contain  the  actual  seismograms.  This  database  format  is 
compatible  with  the  SRO  magnetic-tape  format11  and  the  format  of  the  data  received  from  the 
Datacomputer  through  the  DCREAD  program  (see  Sec.  F below). 

The  general  display  program  will  either  produce  a display  from  an  existing  command  file 
or  create  a default  command  file  for  the  database.  Once  the  initial  display  has  been  made,  the 
program  serves  as  an  Interface  between  the  user  and  the  display,  allowing  the  user  to  type  in  a 
series  of  change  commands.  The  program  will  create  a temporary  display  command  file  from 
these  changes  and  modify  the  display  accordingly. 

The  change  commands  will  allow  the  user  to  rescale  the  data  in  both  x and  y,  move  the 
waveform  left,  right,  up,  or  down,  change  the  separation  of  the  waveforms,  set  up  a cursor, 
align  waveforms  at  the  cursor  location,  and  set  up  a marker  at  the  cursor  location.  All  changes 
will  be  plotted  without  clearing  the  screen,  unless  a clear  command  has  been  given.  The  old 
command  file  will  then  be  saved  in  a temporary  file  and  replaced  by  the  modified  command  file. 
There  will  also  be  a change  command  which  will  allow  the  user  to  return  the  display  to  this 
previous  file. 

The  input  to  the  display  subroutine  will  be  the  display  command  file  which  can  have  an  un- 
limited number  of  lines  of  varying  lengths.  The  commands  are  in  a free  format  where  the  first 
three  items  must  be  the  command,  the  waveform  number  the  command  applies  to,  and  a flag 
which  determines  whether  or  not  the  item  should  be  visible.  The  remaining  parameters  in  each 
line  of  the  command  specify  the  details  of  the  display.  The  display  subroutine  will  evaluate  each 
line  of  this  command  file  and  create  graphics  commands  to  produce  w single  picture.  A large 
portion  of  the  coding  for  this  subroutine  has  been  done.  Figure  IV-li  shows  the  display  that  is 
the  result  of  processing  the  following  command  file: 

waveform  01  1 0.  -1500  1000  2 .5  2000  250 

waveform  02  1 0.  -1500  0 2 .5  2000  250 

waveform  03  1 0.  -1500  -1000  2 .5  2000  250 

station  all  1 -500  200 

Inst  all  1 -500  120 

date  aU  1 -500  40 

time  all  1 -500  -40 

file  all  1 -500  -120 

comment  01  1 -500  -200 

text  none  1 -650  -1500  "TEST  OF  WAVEFORM  DISPLAY  PACKAGE" 
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The  "waveform"  command  defines  where  a waveform  will  appear  on  the  screen  and  sets  up 
the  graphic  commands  to  actually  plot  the  seismogram.  There  must  be  a "waveform"  command 
for  each  waveform  file  that  is  to  be  displayed.  The  remaining  commands  in  the  display  command 
file  r-  quire  an  x and  y offset  from  the  waveform  positions.  The  "station,"  "inst,"  "file,"  and 
"comment"  c .mmands  will  retrieve  the  station  code,  instrument  response  and  component  infor- 
mation, the  gram  file  name,  and  the  comment,  respectively,  from  the  index  file  for  this  database 
and  plot  them  on  the  screen.  If  date  and  time  commands  are  given,  the  display  subroutine  will 
calculate  the  date  and  time  of  the  first  data  point  in  the  display  window  by  using  the  date,  time, 
and  sampling  interval  parameters  from  the  index  file.  Since  the  text  command  is  used  for  la- 
beling and  the  information  is  not  necessarily  associated  with  one  particular  gram,  "none"  is  used 
instead  of  a gram  number. 

M.  F.  Patton 
L.  J.  Turek 


F.  DATACOMPUTER -RELATED  UNIX  SOFTWARE 

A number  of  tools  have  been  developed  to  give  us  easy  access  to  the  seismic  data  on  the 
Datacomputer.  First,  we  now  have  a general  Datacomputer  subroutine  package  (DCS)  which  pro- 
vides the  basic  Datacomputer  functions  to  any  C program.  These  functions  include: 

dcinit  - open  a Datacomputer  connection 
dcdl  - send  datalanguage  to  the  Datacomputer 
dcstart  - start  a data  transfer 

dcopen  --  define  a separate  connection  to  the  Datacomputer 
dc-quit  - close  the  Datacomputer  connection 

Using  the  DCS  package,  we  have  written  a program  (DCREAD)  to  get  SRO  data  from  the 
Datacomputer  and  store  them  in  the  standard  UNIX  waveform  database.  A waveform  database 
consists  of  an  alphanumeric  index  file  which  contains  one  line  describing  each  segment  of  data, 
and  a series  of  binary  files  which  contain  the  actual  seismograms. 

The  format  of  a DCREAD  request  is  the  same  as  the  format  of  the  standard  waveform  index 
file.  It  contains  a series  of  lines,  each  line  describing  a desired  segment  of  data  by  Its  station 
i^de  start  time,  number  of  samples,  SP  or  LP  response,  and  component.  DCREAD  processes 
each  line  individually;  if  a data  request  cannot  be  satisfied,  it  gives  an  error  message  and  moves 
on  to  the  next  line.  LP  requests  will  always  result  in  a single  output  waveform,  but  SP  requests 
will  return  a separate  output  waveform  for  each  data  segment  found  within  the  requested  window. 

Because  DCREAD  does  not  require  any  operator  interaction,  it  is  practical  to  construct 
quite  large  data  requests  and  to  let  the  program  run  an  long  as  necessary  to  retrieve  all  the  data. 

There  are  a few  inconveniences  involved  in  using  DCREAD  for  all  of  our  data  retrieval. 

First,  DCREAD  must  access  the  raw  NLPF  and  NSPF  data  files.  With  the  current  configuration 
of  the  Datacomputer,  these  files  are  on-line  for  only  a very  short  time  (usually  less  than  a month). 
Second,  since  the  data  must  be  specified  by  station  name  and  arrival  time,  the  user  has  the  re- 
sponsibility of  finding  an  event  list  and  computing  the  arrival  times  to  be  used  in  the  DCREAD 
request. 

To  circumvent  these  problems,  we  plan  to  create  our  own  Datacomputer  files  which  will  con- 
tain event  information  from  the  PESF  (Preliminary  Event  Summary  File),  the  computed  SRO 
arrival  times,  and  the  associated  waveform  segments.  These  files  will  stay  on-line  permanently, 
rather  than  being  cycled  off-line  along  with  the  raw  data,  and  the  waveform  data  can  be  accessed 
via  searches  directly  on  the  event  parameters. 
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The  procedure  for  netting  up  these  files  will  include  several  distinct  steps: 

(1)  Retrieve  the  PESF  event  information  for  a given  time  period  from  the 
Datacomputer  to  a local  UNIX  file. 

(2)  Add  to  this  list  any  additional  events  that  we  know  about  from  other 
sources. 

(3)  Retrieve  the  SP  detection  windows  for  the  time  period.  (Start  and  stop 
times  only.) 

(4)  Compute  arrival  windows  for  each  event  at  each  SRO  site  for  both  SP 
and  LP  signals.  (The  SP  window  will  match  the  SP  detection  starting 
nearest  to  the  arrival  time.  The  LP  window  will  start  with  the  origin 
time  and  vary  in  length  depending  on  the  site's  distance  from  the  event.) 

(5)  Move  data  for  those  windows  from  the  raw  data  files  to  a Lincoln 
Waveform  File  (LWF). 

(6)  Add  the  event  information  and  computed  arrival  times  to  the  Lincoln 
Event  Summary  File  (LESF). 

The  program  that  will  carry  out  these  steps  is  being  designed  with  recovery  procedures  to 
Insure  easy  restarting  after  network  failures.  Most  of  the  design  work  for  the  program  is  com- 
plete. and  software  for  steps  (1)  and  (3)  and  part  of  steps  (4)  and  (5)  has  been  written. 

L.  J.  Turek 
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Fig.  IV-i.  Seismometer  broadband 
response  to  acceleration  input, 
(a)  Amplitude  response  relative  to 
DC  response;  (b)  phase  response. 
The  -1  factor  In  transfer  function 
corresponds  to  + ir  radians  phase 
shift. 


Fig.  IV-2.  Seismometer  filtered 
response  to  displacement  input, 
(a)  Amplitude  response  relative 
to  peak  value;  (b)  phase  response. 
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Fig.  IV-3.  Overall  LP  SRO  response 
to  ground  displacement  with  and  with- 
out anti-alias  filter,  (a)  Amplitude 
response;  (b)  phase  response. 


Fig.  IV-4.  Overall  SP  SRO  response 
to  ground  displacement  with  and  with- 
out anti-alias  filter,  (a)  Amplitude 
response;  (b)  phase  response. 


Fif.  IV-9.  First  arrivals  determined  by  computer.  Vertical  dotted  line 
is  computed  arrival  time,  and  number  to  right  of  it  is  estimate  of  standard 
error  of  arrival  in  seconds.  For  each  trace,  1 sec  of  noise  before  and 
i sec  of  signal  after  arrival  time  is  shown.  Tick  marks  are  0.1  sec  apart. 
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Pig.  IV-tO.  Using  a tentative  location  to  constrain  arrival  time. 
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Fig.  SV-li.  Displayed  to  loft  of  each  seismogram  are  the  station  code, 
instrument  cods,  date  and  time  of  the  first  sample,  the  gram  file  name, 
and  comment. 
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ALPA 

ARPANET 

Alaakan  Long -Period  Array 

DARPA  Computer  Network 

DADS 

DCS 

Data  Analyala  and  Diaplay  Syatem 
Datacomputer  Subroutine  Package 

EDR 

Earthquake  Data  Reporta 

IFFT 

ISC 

Inverae  Fast  Fourier  Tranaform 
International  Selamologtcal  Center 

J-B 

Jeffreya-Bullen  (Tablea) 

LASA 

LP 

LRSM 

Large  Aperture  Selamic  Array 

Long  Period 

Long  Range  Selamic  Meaaurementa 

NCER 

NTS 

National  Center  for  Earthquake  Reaearch 

Nevada  Teat  Site 

PDE 

PESF 

Preliminary  Determination  of  Eplcentera 
Preliminary  Event  Summary  File 

SATS 

SDP 

SNR 

SP 

SRO 

Semiannual  Technical  Summary 

Signal  Diaplay  Package 

Signal -to -Noiae  Ratio 

Short  Period 

Selamic  Reaearch  Obaervatory 

USGS 

United  Statea  Geological  Survey 

WWSSN 

World  Wide  Standard  Selamograph  Network 
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